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ABSTRACT 


This manual describes the organization and flow of data and analysis modules in the computer 
code, HITCAN (High Temperature Composite ANalyzer). HTTCAN is a general purpose 
computer program for predicting nonlinear global structural and local stress-strain response of 
arbitrarily oriented, multilayered high temperature metal matrix composite structures. This manual 
describes the architecture of the HITCAN code, followed by the listing of subroutines and calling 
tree, data storage scheme, file system, and a dictionary of code terminology. The primary 
intention of the manual is to familiarize the user with some of the computer program related 
issues so as to facilitate maintenance/modification/updates of the HITCAN computer code. 
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CHAPTER 1 

INTRODUCTION 

HITCAN (High Temperature Composite ANalyzer) is a general purpose computer 
program for predicting global structural and local stress-strain response of high 
temperature metal matrix composite structures. This document describes the 
internal architecture of the code. 

The HITCAN computer program is a combination of three computer programs 
developed in-house at the NASA Lewis Research Center (LeRC) . They are: 

- a dedicated mesh generator; COBSTRAN (Reference 1) 

- a finite element structural analysis code; MHOST (Reference 2) 

- a metal matrix composite analyzer; METCAN (Reference 3). 

HITCAN itself is made up of approximately 112 subroutines consisting of 

16.000 lines of coding. COBSTRAN, METCAN, and MHOST consist of approximately 

7.000 10000, and 51000 lines, respectively. HITCAN is written in FORTRAN 77. 
At the present time HITCAN has been configured and executed on both the CRAY XMP 
and YMP computers at NASA LeRC. 

This document describes: 

- the control structure of the program; 

- the data storage scheme and the memory allocation procedure ; 

- the file management facilities. 

The architecture of the program is described in Chapter 2. The HITCAN calling 
tree and a brief description of each subroutine can be found in Chapter 3. The 
common blocks and the memory allocation procedure is reviewed in Chapter 4. The 
file system is described in Chapter 5. Finally, in Chapter 6, a 
dictionary of most of the variables used in HITCAN is provided. 
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CHAPTER 2 

PROGRAM ARCHITECTURE 

As mentioned previously, HITCAN is the synthesis of three in-house computer 
programs COBSTRAN, MHOST, and METCAN, Figure 2.1, illustrates what the task is 
of each program. As seen in the figure, COBSTRAN determines the finite element 
model, ply lay-up at each node, and the Interpolated nodal loads. MHOST 
determines the structural response of the structure and METCAN calculates the 
laminate properties, ply stresses, and constituent stresses. Figure 2.2 shows 
the how the three programs fit together in HITCAN. As seen in this figure 
COBSTRAN needs to be called only once. MHOST and METCAN, which are part of an 
incremental nonlinear solution scheme, must be called repeatedly. In this 
scheme, METCAN and MHOST continually feed information back and forth, until 
convergence is reached. 

To see the actual connections between COBSTRAN, MHOST, and METCAN, the 
architecture of HITCAN needs to be examined. Figure 2.3 shows the flow chart of 
the executive control module . This module initializes the program FORTRAN units , 
stores the user's input deck for later reading, and initializes the program 
options. The options selected are used to control program flow through the three 
second level executive subroutines HSOLID, HPLATE, and S3DS0L. These secondary 
executive subroutines initialize storage arrays, process the input files, 
generate the finite element model, and perform either an incremental static 
analysis or a dynamic analysis. The general flow diagram shown in Figure 2.4 
corresponds to each one of the three second level executive subroutines . This 
can be seen in their flow charts (See Sections 2.1 to 2.3). Note that the module 
COBSTRAN is directly called from these subroutines. MHOST and METCAN are called 
indirectly from the second level subroutines via the subroutines NLINER and 
DYNMCH. Since the dynamic analysis is not complete, only the interaction between 
MHOST and METCAN in subroutine NLINER will be described. For a summary of 
subroutine NLINER, see Section 2.4. 
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HITCAN: AN INTEGRATED APPROACH FOR HIGH TEMPERATURE 
COMPOSITE STRUCTURAL ANALYSIS 
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The secondary executive subroutine HSOLID controls the generation and the 
analysis of solid structures using plane stress, plane strain, or plate elements. 
This executive subroutine is described in Section 2.1. The function of the 
subroutine HP1ATE is to control the generation and analysis of hollow structures 
using plate elements. A description of this subroutine is given in Section 2.2. 
The subroutine S3DS0L is used to control the generation and analysis of solid 
structures using eight-node solid elements. This subroutine is described in 
Section 2.3. 

A dynamic core allocation strategy has been implemented in the HITCAN 
program. There are three storage arrays used in HITCAN. The array A is used for 
storage in METCAN, the array B is used in HITCAN, and the array IWORK is used as 
the storage array in MHOST. The variables MAXL.EM, NASIZE, and ISIZE are the 
length of these arrays, respectively. These variables control the core memory 
requirements for the HITCAN program. The three arrays A, B, and IWORK and the 
variables MAXLEN, ISIZE, and NASIZE are set in the executive control module. If 
the number of words required to do an analysis exceeds the size of the array B, 
execution will stop and a message will be printed with the number of additional 
words needed. The number of words of allowed storage should be increased and the 
code re-compiled, so that the analysis can continue. 
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NOTE: ( ) Denotes a Subroutine Name. 
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2.1 SUBROUTINE HSOLID 

This subroutine is a second level executive subroutine, which controls the 
generation and the analysis of solid structures using plane strain, plane stress, 
or plate elements. The general flow diagram of this subroutine is shown in 
Figure 2.3. This subroutine performs the following functions : 

1. Initialization of the storage array B. This is rendered through the 
subroutines SCPEM1 , SCPEM2 , INITP3, and SCTEMW. 

2. Processing of user input file. This function is done through the 

subroutines RINPAR, SREAD1, SREAD2, RINPLY, SREAD3, RINFOR, DYNMIN, 
RIBOUN, and OUTPTC. -v>. — • 

3. Generation of the finite element model. This job is carried out in the 
subroutine PREPSO. This subroutine is a modified version of COBSTRAN, 
see Reference 4. 

4. Incremental static or dynamic analysis. These two functions are 
performed in the subroutines NLINER and DYNMCH. A description of the 
subroutine NLINER can be found in Section 2.4. 

A detailed flow chart of this subroutine is shown in Figure 2.5. 
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NOTE: ( ) Denotes a Subroutine Name. 
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2.2 SUBROUTINE HFLATE 

This subroutine is a second level executive program, which controls the 
generation and the analysis of hollow structures using a plate element. This 
subroutine follows the general flow diagram shown in Figure 2.3. The functions 
performed in this subroutine are : 

1. Initialization of the storage array B. This is rendered through the 
subroutines HCPEMl, HCPEM2 , INITP3, and HCTEMW. 

2. Processing of user input file. This function is done through the 
subroutines RINPAR, RINSPA, HREAD1 , HREAD2 , RINPLY, HREAD3 , RINFOR, 
DYNMIN, RIBOUN, and OUTPTC. 

3. Generation of the finite element model. This job is carried out in the 
subroutines PREPSH and PREPSP. These subroutines taken together are a 
modified version of COBSTRAN, see Reference 4. 

4. Incremental static or dynamic analysis. These two functions are 
performed in the subroutines NLINER and DYNMCH. A description of the 
subroutine NLINER can be found in Section 2.4. 

The flow chart of this subroutine is shown in Figure 2.6. 
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NOTE: ( ) Denotes a Subroutine Name. 
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NOTE: ( ) Denotes a Subroutine Name. 
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NOTE: ( ) Denotes a Subroutine Name. 
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2.3 SUBROUTINE S3DSOL 

This subroutine is a second level executive program, which controls the 
generation and the analysis of solid structures using the eight-node solid 
element. This subroutine follows the general flow diagram shown in Figure 2.3. 
The functions performed in this subroutine are: 

1. Initialization of the storage array B. This is rendered through the 
subroutines BCPEM1, BCPEM2, BCPEM3 , and INITP3. 

2. Processing of user input file. This function is done through the 
subroutines RINPAR, BREAD 1 , BREAD 2 , READML, BREAD 3 , RINPLY, RINLAY, 
BREAD4 , RINFOR, DYNMIN, RIBOUN, and OUTPTC. 

3. Generation of the finite element model. This job is carried out in the 
subroutines PCOBST, COBSTRAN, NASINT, and METINT. The subroutine PCOBST 
creates the input file to COBSTRAN, see Reference 4. The subroutines 
NASINT and METINT interpret the files from COBSTRAN. 

. Incremental static or dynamic analysis. These two functions are 
performed in the subroutines NLINER and DYNMCH. A description of the 
subroutine NLINER can be found in Section 2.4. 

The flow chart of this subroutine is shown in Figure 2.7. 
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NOTE: ( ) Denotes a Subroutine Name. 
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NOTE: ( ) Denotes a Subroutine Name. 


Chapter 2 


FLOW CHART OF SUBROUTINE S3DSOL March, 1992 
FIGURE 2.7 (continued) 


ARCH19 




HITCAN Programmer's Manila 1 - Version 1.0 



NOTE: ( ) Denotes a Subroutine Name. 
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2.4 SUBROUTINE NLINER 

This subroutine controls the sequence of execution for the incremental static 
analysis. The scheme used in HITCAN is the incremental direct iteration method. 
The library modules MHOST and METCAN play an integral part in this method. 
Before examining the roles of MHOST and METCAN in the direct iteration method, 
the scheme itself needs to be described. 

The solution scheme in a general form, assuming a single degree-of -freedom 
system, works as follows. 

Initially, the material properties are calculated based on the nodal 
temperatures and the nodal stresses obtained from the previous load 
increment . 

The incremental displacements are calculated from the equation 
AD 1 - [K°]' 1 AF, 

where [K°] is the initial stiffness matrix and AF is the 
incremental load. 

Repetition of this process can be written as 
AD 11 - [K^j^AF. 

Here the stiffness matrix [K 11-1 ] Is calculated based on the material 
properties of the previous iteration. The material properties are 
functions of the nodal temperatures and stresses. 

This process Is terminated when the error becomes sufficiently small, i. 
e. , 

j j AD 11 - AD n_1 ! j/j |AD n_1 - AD 11 ' 2 1 | < tolerance. 

This process is shown graphically in Figure 2.8. 
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A general flow diagram of the subroutine NLINER is shown in Figure 2.9. 
Figure 2.10 shows how MHOST and METCAN are incorporated into the solution 
scheme. At the beginning of each iteration METCAN, block 1, is called. Using 
the constituent properties, the ply stresses, the constituent stresses, and 
the new temperature, METCAN determines the new laminate properties for each 
node. This means METCAN must be called for each node. MHOST (block 2) using 
the laminate properties from METCAN and the nodal loads , determines the 
incremental nodal displacements and the incremental nodal stress resultants . 
The new incremental displacements are then compared with the previous 
incremental displacements to check for convergence. If the solution fails to 
converge, METCAN (block 1) determines new laminate properties. However, if 
COTXVer gence is achieved, the nodal displacements are added to the incremental 
displacements (block 3) and METCAN is then called again to determine the new 
incremental ply stresses and constituent stresses (block 4) . These stresses 
are then added to the cumulative stresses in block 5. It should be pointed 
out that METCAN is called for each node so that the effects at the 
micromechanics level can be captured. 

A detailed flow chart of this subroutine is shown in Figure 2.11. 
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NOTE: ( ) Denotes a Subroutine Name. 
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NOTE: ( ) Denotes a Subroutine Name. 
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CHAPTER 3 

SUBROUTINE LISTING AND CALLINGTREE 

The HITCAN calling tree and brief descriptions of the subroutines included 
in the HITCAN code are given in this chapter. Figure 3.1 shows the HITCAN 
calling tree. 

A brief description of each routine is given below. The subroutine names are 
sorted in alphabetical order. Note that almost all subroutines written for the 
HITCAN computer code are self- documented and further information can be obtained 
directly from the source listing. 

BANKRH 

Reads in the fiber, matrix, and interface properties. Called from the 
subroutines PREPSH, PREPSO, and S3DS0L. 

BCPEM1 

Sets the pointers NB1 , NB2, and NB3. Called from the secondary executive 
subroutine S3DS0L. 

BCEPM2 

Sets the pointer NB4. Called from the secondary executive subroutine S3DS0L. 


BCPEM3 

Sets the pointers NB5 , NB6 NB39 , NB40. Called from the secondary 

executive subroutine S3DS0L. 

BNORML 

Generates the nodal coordinate system and determines the transformation 
matrix from the global to local coordinate systems. Also computes the normals 
to the nodes. Called from the secondary executive subroutine S3DS0L. 
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BREAD 1 

Reads in the variables NXSPC, XH, and LSECT. Called from the secondary 
executive subroutine S3DSOL. 

BREAD2 

Reads in the array MSECT. Called from the secondary executive subroutine 
S3DS0L. 

BREAD 3 

Reads in the array PIN. Called from the secondary executive subroutine 
S3DS0L. 

BREAD4 

Reads in the arrays TISTPS, TEPR, TPGNP, NPRESS, and PREVAL. Called from the 
secondary executive subroutine S3DS0L. 

BUCKSH 

Generates the Input file for the MHOST module (SOLVEH) , calls SOLVEH, and 
writes the eigenvalues and eigenvectors to the output file. Called from the 
subroutine NLINER. 

CHCKMS 

Determines the number of words required in HITCAN and compares this value to 
the size of the array B. Called from the secondary executive subroutines HSOLID, 
HPLATE , and S3DS0L. 

CHECKD 

Compares the new incremental nodal displacements with the previous 
incremental displacements using the norm. Called from the subroutine NLINER. 
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DETNPL 

Determines the maximum nodal thickness in the model and calculates the 
maximum number of plies required to fill that thickness. Called from the 
secondary executive subroutine S3DS0L. 

DETPAR 

Determines several parameters that are rquired in allocating storage. Called 
from the secondary executive subroutine S3DS0L. 

DYNMIN 

Reads in the data for the dynamic analysis. This data consists of the arrays 
NDISPL, DISPIN, DISPI, NVELO, VELOIN, VELOI , NACCEL, ACELIN, ACCELI , NPERD , and 
PERDLD. Called from the secondary executive subroutines HSOLID, HPLATE , and 
S3DS0L. 

DYNMCH 

Performs the dynamic analysis. Called by the secondary executive subroutines 
HPLATE, HSOLID, and S3DS0L. 

FFREADH 

This subroutine is used to support the free-field input format when reading 
program option cards. Called from the subroutine OPTNH. 

E11C0N 

Writes the element connectivities of the plane strain element to the MHOST 
input file. Called by the subroutine MHPREP. 

E75C0ND 

Writes the element connectivities of the plate element to the MHOST input 
file. Called by the subroutine MTPREP. 
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EL3C0N 

Writes the element connectivities of the plane stress element to the MHOST 
input file. Called by the subroutine MHPREP. 

EL7C0N 

Writes the element connectivities of the eight node solid element plane 
strain to the MHOST input file. Called by the subroutine MHPREP. 

EL75CON 

Writes the element connectivities of the plane strain element to the MHOST 
Input file. Called by the subroutine MHPREP. 


EL3PROP 

Writes the D-matrix to the MHOST input file for the plane stress element. 
Called by the subroutine MHPREP. 

EL7PR0P 

Writes the D-matrix to the MHOST input file for the eight node solid element. 
Called by the subroutine MHPREP. 

E11PR0P 

Writes the D-matrix to the MHOST input file for the plane stress element. 
Called by the subroutine MHPREP. 

E75PROF 

Writes the nodal laminate properties to the MHOST input file for the plate 
element. Called by the subroutine MHPREP. 

ELTYPE 

Sets various parameters based on the type of element used in the analysis. 
This subroutine is called by the executive module. 
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GETIME 

A utility subroutine, whose function is to obtain the cumulative cpu time. 


HCPEM1 

Sets the pointers NH1 , NH2 NH32 , NH33 . Called from the secondary 

executive subroutine HPLATE. 

HCPEM2 

Sets the pointers NH34, NH35, NH45 , NH46. Called from the secondary 

executive subroutine HPLATE. 

HCTEMW 

Sets the pointers NH36T, NH37T NH53T, NH54T. Called from the 

secondary executive subroutine HPLATE. 

HDNPLY 

Determines the maximum wall thickness of the shell and calculates the number 
of plies required to fill the wall thickness. Called from the secondary 
executive subroutine HPLATE. 

HINIT 

Sets the FORTRAN unit numbers. Also opens various files. This subroutine 
is called by the executive module. 

HPLATE 

This subroutine is a second level executive program, which controls the 
generation and the analysis of hollow type structures using a plate element. 

HREAD1 

Reads in the array MWSECT. Called from the secondary executive subroutine 


HPLATE. 
Chapter 3 
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HREAD2 

Reads in the array PIN. Called from the secondary executive subroutine 
HPLATE. 

HREAD3 

Reads in the arrays TISTPS , TEPR, NPRESS, NTNOD, TEPLY, and PREVAL. Called 
from the secondary executive subroutine HPLATE. 

HSOLID 

This subroutine is a second level executive program, which controls the 
generation and the analysis of solid type structures, using plane strain, plane 
stress, or plate elements. 

HTMET 

Calls the HITCAN library module METCAN. Called from the subroutine NLINER. 
INITMS 

Initializes the microstresses, microstress rates, microstrains, ply stresses, 
ply strains , and cumulative global displacements . Called from the subroutine 
NLINER. 

INITP3 

Sets the pointers ND42, ND43 ND83, ND84. Called from the secondary 

executive subroutines HPLATE, HSOLID, and S3DS0L. 

INITSR 

Initializes nodal stress resultants. Called from the subroutine NLINER. 
LOGOH 

Writes to the output file the HITCAN logo. This subroutine is called by the 
executive module . 
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MHPREP 

Writes the input file KMHOST into the module MHOST. Called from the 
subroutine NLINER. 

MIDWAL 

Moves the input points from the shell surface to the shell wall mid- thickness 
line. This subroutine is used in the HPLATE model option. Called by the 
subroutine WPROFIL. 

MODESH 

Generates the input file for the MHOST module (SOLVEH) , calls SOLVEH, and 
writes the frequencies, eigenvalues, and eigenvectors to the output file KOUT. 
Called from the subroutine NLINER. 

MPROP 

Reads the microproperties calculated in METCAN from the post files KP0ST2, 
KP0ST3 , KP0ST4 , and KP0ST5 . Sets the current material flags and writes the 
microproperties to the output file KOUT. Called from the subroutine NLINER. 

MSTRES 

Reads the ply stresses and microstresses calculated in METCAN from the post 
files KP0ST6 , KP0ST7 , KP0ST8 , and KP0ST9 and adds these stresses to the 
cumulative stresses. Also, writes the microstresses to the output file KOUT. 
Called from the subroutine NLINER. 

MTPREP 

Writes the input file KMHOST into the module MHOST. Called from the 
subroutine DYNMCH. 
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NASINT 

Interprets a NASTRAN bulk data deck (generated by COBSTRAN) . This bulk data 
deck contains the finite element model. Called from the secondary executive 
subroutine S3DS0L. 


NLINER 

Performs the incremental static analysis. Called by the secondary executive 
subroutines HPLATE, HSOLID, and S3DS0L. 


NULIN 

A utility routine that zeros out an integer matrix. 

NULREA 

A utility routine that zeros out a real matrix. 

OMODEL 

Writes to the output file KOUT the element connectivities, nodal coordinates, 
and nodal temperatures. Called from the secondary executive module S3DS0L. 


OPTNH 

This subroutine reads the program option cards and sets all logical controls 
in accordance with user choices. Called from the executive control module. 


OTPUTH 


Writes to the output file KOUT, the ply lay-up at each 
element model, load data, and nodal radii of curvature, 
executive control module . 


node, the finite 
Called from the 
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0TJTH10 

Writes to the output file KOUT, the program option cards chosen. Called from 
the executive control module. 

OUT PTC 

Reads in the arrays NPRT, NPRTS, NPRTP , NPPLY, TIMEPN, TIMEMS , and TIMEMB. 
These arrays control the output and also determine when a PATRAN results file is 
to be created, a modal analysis is performed, and buckling analysis is done. 
Called by the secondary executive subroutines HPLATE, HSOLID, and S3DS0L. 

PCOBST 

Generates the input file KCBST, for COBSTRAN. Called from the secondary 
executive subroutine S3DS0L. 

PLYH 

Calculates the ply properties PLH(l) , PLH(2) , PLH(13), PLH(14) . Called 

from the subroutine NLINER. 

PLPROP 

Reads the ply properties from the METCAN post file KP0ST5. Also, determines 
the compliance matrix and the inverse of the compliance matrix. Called from the 
subroutine NLINER. 

PLYTEM 

Sets the ply temperatures. Called from the subroutines PREPSH, PREPSO, and 
PREPSP . 

PMETCN 

Generates the input files KMETCN and KDBANK for the module METCAN. Called 
from the subroutine NLINER. 
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PREPSH 

Cobstran executive control module* Called from tbe secondary executive 
subroutine HPLATE. 

PREPSO 

Cobstran executive control module. Called from the secondary executive 
subroutine HSOLID. 

PREPSP 

Generates that finite element associated with the spars. Called from the 
secondary executive subroutine HPLATE. 

PROPH 

Determines the axial, bending, and coupled stiffness matrices. Called from 
the subroutine NLINER. 

PRTEPS 

Writes the microstrains to the output file KOUT. This subroutine is not used 
at the present time. 

PRTMPR 

Writes the current constituent properties to the output file KOUT. Called 
from the subroutine MPROP. 

PRTMS 

Writes the microstresses to the output file KOUT. Called from the subroutine 
MSTRES . 


Writes the microstress rates to the output file KOUT. Not used at the 
present time. 
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PSTRES 

Write the ply stresses, global displacements, and global reactions to the 
output file KOUT. Called from the subroutine NLINER. 

PSTRSS 

Write the ply stresses, global displacements, and global reactions to the 
output file KOUT. Called from the subroutine DYNMCH. 

READML 

Reads in the nodal coordinates and element connectivities, the arrays KEI and 
GNP, from the input file KREAD. Also, determines the nodes above and below each 
node. Called from the secondary executive subroutine S3DS0L. 

REFLECH 

Temporarily store the user input file in the file KREFL and print the user 
input file, if requested by the user. Called from the executive control module. 

RESTART 

Read/write the restart file KSTART . Called from the subroutine NLINER. 
RIBOUN 

Reads in the boundary conditions and the nodal transformation data from the 
input file KREAD, this includes the arrays NTRANS , TRANG, and NBOUND . Called by 
the secondary executive subroutines HPLATE, HSOLID, and S3DS0L. 

RINFOR 

Reads in the nodal forces and centrifugal loading data from the input file 
KREAD, i. e. , the variables GRIDP1, GRIDP2 , ANGVEL, NCFNOD , NCFDIR, and NCFVAL. 
Called by the secondary executive subroutines HPLATE, HSOLID, and S3DS0L. 
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RINIAY 

Reads in the ply lay-up at each node. This includes the arrays NPNOD and 
MPLY. Called by the secondary executive subroutine S3DS0L. 

RINPAR 

Reads in the parameter data from the input file KREAD. This data consists of 
the variables NDES, NSECT, IGRD, IU, JU, MESH, XBEG, XEND, NXSPAR, NYSPC, MAXNP, 
NETOT, LSECT , NIPL, NOSC, NEYY, NETT, MAXPLY, LMAX, NTISTP, NMECHC , NTHERC, LINC, 
MSTART, MITER, TOL, NEIGV, NSUBD , INCREG, MHITER, RESID, DAMPMS , 

DAMPST, NPLSET, INCDYN , NCFOR , NPRES , NTEMP, NOPERD, NBC, NTR, NDIS, NVEL, NACC, 
and PINTER. Called by the secondary executive subroutines HPLATE, HSOLID, and 
S3DSOL. 

RINPLY 

Reads in ply data, i. e., the arrays PERT, CODEH, MPLY, and NPLY from the 
input file KREAD. Called by the secondary executive subroutines HPLATE, HSOLID, 
and S3DSOL. 

RINSPA 

Reads in the spar description data from the input file KREAD. This data 
consists of the arrays NSPDES, NXSPC, XH, and SY. Called by the secondary 
executive subroutine HPLATE. 

S3DS0L 

This subroutine is a second level executive program, which controls the 
generation and the analysis of solid type structures, using the eight node solid 
element . 

SCPEM1 

Sets the pointers NM1 and NM2 . Called from the secondary executive 
subroutine HSOLID. 
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SCPEM2 

Sets the pointers NM3, NM4, NM60, NM61. Called from the secondary 

executive subroutine HSOLID. 

SCTEMW 

Sets the pointers NMT61, NMT62 , NMT72, NMT73 . Called from the 

secondary executive subroutine HSOLID. 

SERCH 

Determines the maximum nodal displacement in the PATRAN results file KNODAL. 
Called from the subroutine WPOST. 

SHELLT 

Determines the transformation matrix from the local nodal coordinate system 
to the global coordinate system. This subroutine is used only for the plate 
element. Called from the secondary executive subroutines HPLATE and HSOLID. 

SOLVEH 

Calls the HITCAN library module MHOST. Called from the subroutine NLINER. 
SREAD1 

Reads in the array MSECT from the input file KREAD. Called from the 
secondary executive subroutine HSOLID. 

SREAD2 

Reads in the arrays NXDIV and PIN from the input file KREAD. Called from the 
secondary executive subroutine HSOLID. 
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SREAD3 

Reads in the arrays TISTPS, TEPR, NTNOD, TEPLY, NPRES , and PREVAL from the 
input file KREAD. Called from the secondary executive subroutine HSOLID. 

TANLOC 

Establishes a set of three orthigonal unit vectors VI, V2, and V3 , such that 
V3 lies along the given vector X and V2 is in the plane defined by the vectors 
X and Y, Called from the subroutine NLINER . 

TBGRDH 

Calculates the nodal coordinates, nodal temperatures, nodal pressure loading, 
and the wall thickness at each node from the user input. This subroutine is used 
in the HSOLID model option. Called from the subroutine PREPSO. This subroutine 
is a modified version of the subroutine TBGRD, which is found in the library 
module COBSTRAN (see Reference 1). 

TBG01H 

Divides a solid structure along the X-axis and Y-axis into the appropriate 
divisions. Called from the subroutine TBGRDH. This subroutine is a modified 
version of the subroutine TBG01, which is found in the library module COBSTRAN 
(see Reference 1). 

TBPLY 

Determines the number of plies required at each node for the HSOLID model 
option. Called from the subroutine PREPSO. This subroutine is a modified 
version of the subroutine TBGRD, which is found in the library module COBSTRAN 
(see Reference 1) . 

WDIS 

Writes the global displacements to the PATRAN nodal results file KNODAL. 
Called from the subroutine WPOST. 
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WINDH 

Generates the element connectivities for the HPLATE model option. Called 
from the subroutine PREPSH. This subroutine is a modified version of the 
subroutine WIND, which is found in the library module COBSTRAN (see Reference 1) . 

WNOD 

Writes the stress resultants to the PATRAN nodal results file KNODAL. Called 
from the subroutine WPOST. This subroutine is not used at the present time. 

WOUTH 

Writes to the output file KOUT the element connectivities, nodal coordinates, 
etc. for the HPLATE model option. Called from the subroutines PREPSH and PREPSP. 
This subroutine is a modified version of the subroutine WINOUT, which is found 
in the library module COBSTRAN (see Reference 1) . 

WPLYH 

Determines the ply lay-up at each node for the HPLATE model option. Called 
from the subroutines PREPSH and PREPSP. This subroutine is a modified version 
of the subroutine WINOUT, which is found in the library module COBSTRAN (see 
Reference 1) . 

WPOST 

Controls the generation of the PATRAN results files KNODAL and KTRAN. Called 
from the subroutine NLINER. 

WPUTH 

Determines the nodal coordinates, nodal temperatures, and nodal pressures for 
the HPLATE model option. Called from the subroutine PREPSH. This subroutine is 
a modified version of the subroutine WINPUT, which is found in the library module 
COBSTRAN (see Reference 1). 
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WNEPT 

Writes a PATRAN neutral file. Called by the secondary executive subroutines 
HPLATE , HSOLID , and S3DSOL. 

WPRQFIL 

Moves the grid points from the surface of the structure to the mid-plane of 
the shell. This subroutine is used in the HPLATE model option. Called from the 
subroutine PREPSH. This subroutine is a modified version of the subroutine 
MIDWAL, which is found in the library module COBSTRAN (see Reference 1) . 

WSPARH 

Determines the element connectivities for the spars. Also, finds the 
corresponding shell node for each spar node. Called from the subroutine PREPSP. 
This subroutine is a modified version of the subroutine WSPAR, which is found in 
the library module COBSTRAN (see Reference 1) . 
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CHAPTER 4 

DATA STORAGE SCHEME 


A dynamic core allocation strategy has been implemented in the HITCAN 
program. The array B is used as the working storage. The length of this 
array is the variable NASIZE. The array B is used to store only the 
information used directly in HITCAN. METCAN and MHOST have their own dynamic 
storage areas. In METCAN the array A is used as the working storage, MAXLEN 
is the length of this array. IWORK is the working storage array for MHOST and 
ISIZE is its length. The three arrays A, B, and IWORK and the variables 
MAXLEN, ISIZE, and NASIZE are set in the main program of HITCAN. If the 
number of words required for an analysis exceeds the size of the array B, 
execution will stop and a message will be printed with the number of 
additional words needed. For a description of the dynamic core allocation 
scheme employed in MHOST, the programmer is referred to Reference 5. 

In HITCAN the core allocation is carried out in the subroutines HCPEM1, 
HCPEM2 , HCTEMW, SCPEM1 , SCPEM2 , SCTEMW, BCPEM1 , BCPEM2 , BCPEM3 , and INITP3 . 

The first three of these subroutines are used for the HPLATE model option, the 
next three subroutines are used in the HSOLID model option, the next three are 
used in the S3DS0LID model option and the READ IN THE MODEL option. The last 
subroutine INITP3 is used by each of the model options. 

In Sections 4.1 through 4.3, the pointers used in the dynamic allocation 
and the common blocks they reside in are listed and described. In Section 4.4 
the common blocks used in HITCAN, which do not contain pointers are also 
listed and described. 
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4.1 POINTERS USED IN THE H PLATE MODEL OPTION 

For the HPLA.TE model option the variables used as pointers and the common 
blocks in which they are located are listed below. The following variables 
are initialized in the subroutine HCPEM1 and are located in the common block 
POINT . 


COMMON /POINT/ 

1 

2 

3 

4 

5 

6 

7 

8 


NH1 , NH2 , 
NH8 , NH9 , 
NH15, NH16, 
NH22 , NH23, 
NH29 , NH30 , 
NH36 , NH37 , 
NH43 , NH44 , 
NH50 , NH51 , 
NH57, NH58 , 


NH3 , 
NH10 , 
NH17, 
NH24 , 
NH31 , 
NH38, 
NH45 , 
NH52 , 
NH59 , 


NH4, NH5 , 
NH11 , NH12 , 
NH18 , NH19, 
NH25 , NH26 , 
NH32 , NH33 , 
NH39 , NH40, 
NH46 , NH47 , 
NH53 , NH54 , 
NH60 , NH61, 


NH6 , NH7 , 
NH13, NH14 , 
NH20, NH21, 
NH27 , NH28 , 
NH34, NH35 , 
NH41 , NH42 , 
NH48 , NH49 , 
NH55 , NH56 , 
NH62 , NH63 


NH1 


Points to the array NSPDES storing the spar descriptions , 


NH2 

NH3 


Points to the array NXSPC storing the number of elements between 
output sections . 

Points to array X, which contains the x- coordinate of the output cross 
sections . 


NH4 Points to the array storing spar information, SY. 

NH5 Pointer to the array (PERT) containing ply information. 

NH6 Pointer to the ply data array CODEH. 

NH7 Points to the array MPLY, which stores the ply lay-up for one half the 

thickness of a symmetrical ply lay-up. 

NH8 Points to the array TISTPS, which contains the time at each time step. 

NH9 Points to the array storing the boundary conditions, NBOUND. 

NH10 Pointer to the array containing the initial fiber properties (FPI) . 
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NH11 Pointer to the array storing the initial matrix properties (MPI) . 

NH12 Pointer to the array storing the initial interface properties (DPI) . 

NH13 Pointer to the array storing the initial fiber stresses (FSIGI) . 

NH14 Pointer to the array containing the initial matrix stresses (MSIGI). 

NH15 Pointer to the array containing the initial interface stresses 

(DSIGI) . 

NH16 Pointer to the array storing the initial fiber stress rates (FDOTI) . 

NH17 Pointer to the array storing the initial matrix stress rates (MDOTI) . 

NH18 Pointer to the array storing the initial interface stress rates 

(DDOTI) . 

NH19 Pointer to the array containing the fiber exponents (FTVCI) . 

NH20 Pointer to the array storing the matrix exponents (MTVCI) . 

NH21 Pointer to the array storing the interface exponents (DTVCI) . 

NH22 Points to the array LSP, which stores the spar location on the top and 
bottom surfaces. 

NH23 Points to the array NF, which stores the number of fibers per bundle 
for each ply. 

NH24 Points to the scratch array NT. 

NH25 Points to the array (DF) storing the fiber diameter for each different 
ply. 

NH26 Points to the scratch array NPD. 

NH27 Points to the array containing pressure data, NPRESS. 

NH28 Points to the array ANGVEL, which contains the angular velocity for 

each time step . 

NH29 Points to the array NCFNOD, which contains concentrated force data. 

NH30 Points to the array NCFDIR, which contains concentrated force data. 

NH31 Points to the array CFVAL, which contains concentrated force data. 
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NH32 Not used. 

NH33 Points to the array MWSCT, which contains the number of input points 
for each cross section. 

The following variables are initialized in the subroutine HCPEM2 and are 
located in the common block POINT . 


NH34 

NH35 

NH36 

NH37 

NH38 

NH39 

NH40 

NH41 

NH42 

NH43 

NH44 

NH45 

NH46 


Points to the array (PIN) storing surface geometry. 

Points to the scratch array NPCE. 

Points to the array (TEPR) storing temperature and pressure data. 

Points to the array PREVAL, which contains edge load intensity for 
each time step. 

Points to the scratch array NCORDL. 

Points to the scratch array RAD. 

Points to the array NTRANS , which stores data for nodal coordinate 
transformation. 

Pointer to the array TRANG, which stores data for nodal coordinate 
transformation. 

Points to the array (DISPI) containing the initial nodal 
displacements . 

Points to the array (VELOI) containing the initial nodal velocities. 

Points to the array (ACCELI) containing the initial nodal 
accelerations . 

Points to the array (NPERD) containing periodic load data. 

Points to the array ( PERDU) ) containing periodic load data. 
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The following variables are initialized in the subroutine HCTEMW and are 
located in the common block POINTT . 


COMMON 

/POINTT/ NH23T, 

NH24T, 

NH25T, 

NH26T, 

i 

NH27T, 

NH28T, 

NH29T, 

NH30T, 

2 

NH31T, 

NH32T, 

NH33T, 

NH34T, 

3 

NH35T, 

NH36T, 

NH37T, 

NH38T, 

4 

NH39T, 

NH40T, 

NH41T, 

NH42T , 

5 

NH43T, 

NH44T, 

NH45T, 

NH46T , 

6 

NH47T, 

NH48T, 

NH49T, 

NH50T , 

7 

NH51T , 

NH52T, 

NH53T, 

NH54T 


NH36T Points to the scratch array TTAB. 

NH37T Points to the scratch array TPTT. 

NH38T Points to the scratch array XINT. 

NH39T Points to the scratch array ZINT. 

NH40T Points to the scratch array VDOT. 

NH41T Points to the array ANG containing the ply orientation angle. 

NH42T Points to the scratch array XR. 

NH43T Points to the scratch array YS. 

NH44T Points to the scratch array YINT. 

NH45T Points to the scratch array TPTTAB. 

NH46T Points to the scratch array TP. 

NH47T Points to the scratch array TPPP. 

NH48T Points to the scratch array TJJ . 

NH49T Points to the scratch array TIN. 

NH50T Points to the scratch array TPTIN. 
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NH51T Pointer to the array NTNOD. This array contains temperature data 
information. 

NH52T Pointer to the array TEPLY. This array contains temperature data 
information. 

NH53T Points to the scratch array TINT. 

NH54T Points to the scratch array PPTS. 


All of the model options use the pointers listed below. These variables 
are initialized in the subroutine INITP3 and are located in the common block 
IPOINT. 


COMMON 

/IPOINT/ ND1 , 

ND2, 

ND3 , 

ND4, 

ND5 , 

ND6, 

ND7 , 

1 

ND8 , 

ND9 , 

ND10 , 

NDll , 

ND12 , 

ND13 , 

ND14 , 

2 

ND15 , 

ND16 , 

ND17 , 

ND18 , 

ND19 , 

ND20, 

ND21, 

3 

ND22 , 

ND23 , 

ND24 , 

ND25 , 

ND26 , 

ND27, 

ND28 , 

4 

ND29 , 

ND30, 

ND31 , 

ND32, 

ND33, 

ND34, 

ND35 , 

5 

ND36, 

ND37 , 

ND38 , 

ND39 , 

ND40, 

ND41 , 

ND42 , 

6 

ND43 , 

ND44 , 

ND45 , 

ND46 , 

ND47 , 

ND48 , 

ND49 , 

7 

ND50, 

ND51 , 

ND52 , 

ND53 , 

ND54, 

ND55 , 

ND56 , 

8 

ND57 , 

ND58 , 

ND59 , 

ND60, 

ND61, 

ND62, 

ND63 , 

9 

ND64, 

ND65 , 

ND66 , 

ND67 , 

ND68, 

ND69, 

ND70 , 

+ 

ND71 , 

ND72 , 

ND73 , 

ND74, 

ND75 , 

ND76 , 

ND77 , 

1 

ND78 , 

ND79 , 

ND80, 

ND81 , 

ND82 , 

ND83 , 

ND84 


ND42 Pointer to the array 
duplicate nodes . 


ND43 

Pointer 

to 

the 

array 

ND44 

Pointer 

to 

the 

array 

ND45 

Pointer 

to 

the 

array 


pressures . 



ND46 

Pointer 

to 

the 

array 


(KAD) storing the shell node numbers of the 

(KEI) storing the element connectivities. 
(GNP) storing the nodal coordinates. 

(TPGNP) storing the nodal temperatures and 

(THHF) storing the nodal half thicknesses. 
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Pointer to the array (NLEP) containing the number of plies at each 
node . 

ND48 Pointer to the array (MDES) storing the ply lay-up. 

ND49 Pointer to the array (THERMO) storing the ply temperatures. 

ND50 Pointer to the axial stiffness matrix (AS). 

ND51 Pointer to the coupling stiffness matrix (CS) . 

ND52 Pointer to the bending stiffness matrix (BS) . 

ND53 Pointer to the coefficients of thermal expansion. 

ND54 Pointer to array storing the ply density (DEN) . 

ND55 Pointer to the array TLAY. This array contains the incremental nodal 

temperatures . 

ND56 Pointer to the array PDIFI , which contains the incremental nodal 
pressures on the lower surface. 

ND57 Pointer to the array PDIFO, which contains the incremental nodal 
pressures on the upper surface. 

ND58 Pointer to the total global displacements (DISP) . 

ND59 Pointer to the incremental global displacements (DISVH) . 

ND60 Pointer to the incremental stress resultants (MSH) . 

ND61 Pointer to the incremental stress resultants (NSH) . 

ND62 Pointer to the incremental stress resultants (QSH) . 

ND63 Pointer to the incremental nodal reactions (REACF) . 

ND64 Points to the incremental nodal loads COFORC. 

ND65 Pointer to the array (PLH) containing the ply properties. 

ND66 Pointer to the array (TOLD) containing the previous ply temperatures. 

ND67 Pointer to the scratch array TEMP. 

ND68 Pointer to the array (TNEW) containing the current ply temperatures. 
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ND69 Pointer to the array containing the eigenvalues (EIGNVA) . 

ND70 Pointer to the array storing the eigenvectors (EIGNVE) . 

ND71 Pointer to the scratch array TEMPER. 

ND72 Pointer to the array (PC), which contains the equivalent nodal 

composite properties. 

ND73 Pointer to the array STRMOM. This array either contains the nodal 
stresses or stress resultants. 

ND74 Pointer to the array PRESUR. This array contains the incremental 
pressure loading. 

ND75 Points to the array (ETRAN) . This array contains nodal transformation 
matrix from global to local coordinates. 

ND76 Pointer to the array (THRSTN) , which contains the thermal strain 
calculated in METCAN. 

ND77 Pointer to the array DISPD. This array contains the initial 
displacements . 

ND78 Pointer to the array VELOD. This array contains the initial 
velocities . 

ND79 Pointer to the array ACCELD. This array contains the initial 
accelerations . 

ND80 Points to the array TOTFOR. 

ND81 Points to the array TOTRED. 

4.2 POINTERS USED IN THE HSOLID MODEL OPTION 


For the HSOLID model option the variables used as pointers and the common 
blocks in which they are located are listed below. The following variables 
are initialized in the subroutine SCPEM1 and are located in the common block 
POINTS . 
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COMMON 

/POINTS/ NM1 , 

NM2 , 

NM3 , 

NM4, 

NM5 , NM6 , 

NM7 , 

1 

NM8 , 

NM9 , 

NM10 , 

NM11, 

NM12, NM13 , 

NM14, 

2 

NM15 , 

NM16 , 

NM17 , 

NM18 , 

NM19, NM20 , 

NM21 , 

3 

NM22 , 

NM23 , 

NM24, 

NM25 , 

NM26 , NM27 , 

NM28 , 

4 

NM29 , 

NM30 , 

NM31 , 

NM32 , 

NM33 , NM34, 

NM35 , 

5 

NM36 , 

NM37 , 

NM38 , 

NM39 , 

NM40 , NM41 , 

NM42 , 

6 

NM43 , 

NM44 , 

NM45 , 

NM46 , 

NM47 , NM48 , 

NM49 , 

7 

NM50, 

NM51 , 

NM52 , 

NM53 , 

NM54, NM55 , 

NM56 , 

8 

NM57 , 

NM58 , 

NM59 , 

NM60 , 

NM61 , NM62 , 

NM63 


NM1 Pointer to the array (MSECT) storing the number of input points for 
each cross-section. 

NM2 Pointer to the array PIN. This array contains the coordinates of the 
input points, which describe the surface geometry. 


The following variables are initialized in the subroutine SCPEM2 and are 
located in the common block POINTS . 

NM3 Pointer to the scratch array NCOOR. 

NM4 Pointer to the array (PERT) containing ply information. 

NM5 Pointer to the ply data array CODEH. 

NM6 Pointer to the array NXDIV, which contains the number of elements 

between output sections. 

NM7 Pointer to the scratch array THWAL. 

NM8 Pointer to the array MPLY. This array stores one half of the ply 

order for a symmetric ply lay-up. 

NM9 Pointer to the array NPLY. This array stores one half of the ply 
order for an unsymmetrical ply lay-up. 

NM10 Pointer to the scratch array NPCE. 

NM11 Pointer to the array containing the initial fiber properties (FPI). 

NM12 Pointer to the array storing the initial matrix properties (MPI) . 
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NM13 Pointer to the array storing the initial interface properties (DPI) . 

NM14 Pointer to the array storing the initial fiber stresses (FSIGI) . 

NM15 Pointer to the array containing the initial matrix stresses (MSIGI) . 

NM16 Pointer to the array containing the initial interface stresses 

(DSIGI) . 

NM17 Pointer to the array storing the initial fiber stress rates (FDOTI) . 

NM18 Pointer to the array storing the initial matrix stress rates (MDOTI) . 

NM19 Pointer to the array storing the initial interface stress rates 

(DDOTI) . 

NM20 Pointer to the array containing the fiber exponents (FTVCI) . 

NM21 Pointer to the array storing the matrix exponents (MTVCI) . 

NM22 Pointer to the array storing the interface exponents (DTVCI) . 

NM23 Pointer to the array (NF) containing the number of fibers per bundle 
for each ply. 

NM24 Pointer to the scratch array NT. 

NM25 Pointer to the array (DF) , which stores the fiber diameters. 

NM26 Pointer to the scratch array NPD . 

NM27 Not used. 

NM28 Pointer to the array storing the time at each time step (TISTPS) . 
NM29 Pointer to the rotational speed (ANGVEL) . 

NM30 Pointer to the location of the concentrated forces (NCFNOD) . 

NM31 Pointer to the array containing the direction of the concentrated 
forces (NCFDIR) . 

NM32 Pointer to the value of the concentrated forces (CFVAL) . 

NM33 Not used. 
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NM34 Pointer to the array NBOUND. This array stores the boundary 
conditions . 

NM35 Pointer to the array storing edge load data (NPRESS). 

NM36 Pointer to the array PREVAL. This array stores the intensity of the 

edge loading at each time step. 

NM37 Pointer to the array storing nodal transformation data (NTRANS) . 

NM38 Pointer to the array TRANG. This array stores the angle of rotation. 

NM39 Pointer to the array storing the initial displacements (DISPI) for 

direct time integration. 

NM40 Pointer to the array storing the initial velocities (VELOI) for direct 
time integration. 

NM41 Pointer to the array storing the initial accelerations (ACCELI) for 
direct time integration. 

NM42 Pointer to the array (NPERD) containing periodic load data. 

NM43 Pointer to the array (PERDLD) containing periodic load data. 

NM55 Pointer to the scratch array NCORDL. 

NM56 Pointer to the scratch array RAD. 

NM57 Pointer to the scratch array XINT, 

NM58 Pointer to the scratch array XOUT. 

NM59 Pointer to the scratch array YINT. 

NM60 Pointer to the scratch array RX. 

NM61 Pointer to the scratch array KDES. 

The following variables are initialized in the subroutine HCTEMW and are 
located in the common block POINTT . 
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COMMON /PINTTS/ 

1 

2 

3 


NMT61, NMT62 , 
NMT65 , NMT66, 
NMT69 , NMT70, 
NMT73 , NMT74, 


NMT63 , NMT64 , 
NMT67 , NMT68 , 
NMT71 , NMT72 , 
NMT75 , NMT76 


NMT61 Points to the scratch array TEPR. 

NMT62 Points to the scratch array TIN. 

NMT63 Points to the scratch array TP. 

NMT64 Points to the scratch array TJ . 

NMT65 Points to the scratch array TT. 

NMT66 Points to the scratch array VDOT. 

NMT67 Points to the scratch array V2D0T . 

NMT68 Points to the scratch array YOUT. 

NMT69 Points to the scratch array RY. 

NMT70 Points to the array NTNOD . This array contains temperature data 
information. 

NMT71 Points to the array TEPLY. This array contains temperature data 
information. 

NMT72 Points to the scratch array TINT. 

NMT73 Pointer to the scratch array PPTS . 


4.3 POINTERS USED IN THE S3DSOLID MODEL OPTION 


For the S3DS0LID model option and the READ IN THE MODEL option the 
variables used as pointers and the common blocks in which they are located are 
listed below. The following variables are initialized in the subroutine 
BCPEM1 and are located in the common block POINTS. 
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COMMON 

/POINTB/ NB1 , 

NB2 , 

NB3 , 

NB4, 

NB5 , 


1 

NB6 , 

NB7 , 

NB8 , 

NB9 , 

NB10 , 


2 

NB11 , 

NB12 , 

NB13 , 

NB14, 

NB15 , 


3 

NB16 , 

NB17 , 

NB18 , 

NB19 , 

NB20 , 


4 

NB21 , 

NB22 , 

NB23 , 

NB24, 

NB25 , 


5 

NB26 , 

NB27 , 

NB28 , 

NB29 , 

NB30 , 


6 

NB31 , 

NB32 , 

NB33 , 

NB34, 

NB35 , 


7 

NB36 , 

NB37 , 

NB38 , 

NB39 , 

NB40 , 


8 

NB41 , 

NB42 , 

NB43 , 

NB44, 

NB45 , 


9 

NB46 , 

NB47 , 

NB48 , 

NB49, 

NB50 , 


+ 

NB51 , 

NB52 , 

NB53 , 

NB54 , 

NB55 


NB1 

Pointer to the array 
output sections. 

(NXSPC) 

storing 

the number 

of elements 

between 

NB2 

Pointer to the array 
output section. 

PIN. This array contains 

the x- ordinate of each 

NB3 

Pointer to the array 

(LSECT) 

storing 

the number 

of sections 

for each 


input plane . 

NB4 Points to the array MSECT. This array stores the number of input 
points for each section. 

NB5 Pointer to the array PIN. This array contains the coordinates of the 
input points, which describe the surface geometry. 

NB6 Points to the array (TEPR) containing temperature data. 

NB7 Pointer to the array (PERT) containing ply information. 

NB8 Pointer to the ply data array CODEH. 

NB9 Pointer to the array MPLY. This array stores one half of the ply 

order for a symmetric ply lay-up. 

NB10 Pointer to the array NPLY. This array stores one half of the ply 
order for an unsymmetrical ply lay-up. 

NB11 Pointer to the array storing the time at each time step (TISTPS) . 

NB12 Pointer to the array NBOUND . This array stores the boundary 

conditions . 

NB13 Pointer to the array containing the initial fiber properties (FPI) . 
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NB14 Pointer to the array storing the initial matrix properties (MPI) . 

NB15 Pointer to the array storing the initial interface properties (DPI) . 

NB16 Pointer to the array storing the initial fiber stresses (FSIGI) . 

NB17 Pointer to the array containing the initial matrix stresses (MSIGI) . 

NB18 Pointer to the array containing the initial interface stresses 

(DSIGI) . 

NB19 Pointer to the array storing the initial fiber stress rates (FDOTI). 

NB20 Pointer to the array storing the initial matrix stress rates (MDOTI) . 

NB21 Pointer to the array storing the initial interface stress rates 

(DDOTI) . 

NB22 Pointer to the array containing the fiber exponents (FTVCI) . 

NB23 Pointer to the array storing the matrix exponents (MTVCI) . 

NB24 Pointer to the array storing the Interface exponents (DTVCI). 

NB25 Pointer to the array (NF) containing the number of fibers per bundle 
for each ply. 

NB26 Pointer to the array (DF) , which stores the fiber diameters. 

NB27 Pointer to the scratch array NPD . 

NB28 Pointer to the rotational speed (ANGVEL) . 

NB29 Pointer to the location of the concentrated forces (NCFNOD) . 

NB30 Pointer to the array containing the direction of the concentrated 

forces (NCFDIR) . 

NB31 Pointer to the value of the concentrated forces (CFVAL) . 

NB32 Pointer to the array storing edge load data (NPRESS). 

NB33 Pointer to the array PREVAL. This array stores the intensity of the 

edge loading at each time step. 

NB34 Pointer to the array storing nodal transformation data (NTRANS) . 
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NB35 Pointer to the array TRANG. This array stores the angle of rotation. 

NB36 Pointer to the array storing the initial displacements (DISPI) for 

direct time integration. 

NB37 Pointer to the array storing the initial velocities (VELOI) for direct 
time integration. 

NB38 Pointer to the array storing the initial accelerations (ACCELI) for 
direct time integration. 

NB39 Pointer to the array (NPERD) containing periodic load data. 

NB40 Pointer to the array (PERDLD) containing periodic load data. 


4.4 COMMON BLOCKS USED IN HITCAN 


The following additional common blocks are used in HITCAN. After a 
listing of each named common block the variables contained in the common block 
are described, where appropriate. Note that if a variable in a common block 
is no longer used in the program this variable is not defined. 


The common blocks HUNITS and CMETCN contain the variable names of the 
Fortran units used in HITCAN. These variables and their assigned unit numbers 


are described in Chapter 6 . 


COMMON 

/HUNITS/ 

KDES1 , 

KKEI , 

1 


KCBST , 

KCBMT, 

2 


KTRAN , 

KPOST , 

3 


KPROPI , 

KSTR, 

4 


KDISP , 

KMSTR, 

5 


KMHOST, 

KSTART , 

COMMON 

/CMETCN/ 

KPOSTO , 

KP0ST1 , 

1 


KPOST5 , 

KP0ST6 , 

2 


KMETCN 



Chapter 4 


KBULK, 

KREAD, 

KOUT, 

KPROP , 

KCBFE , 

KDBANK, 

KMHFOR, 

KREFL, 

KPLY, 

KLAM, 

KSIG, 

KNODAL, 

KTRANF, 

KCBOUT , 

KTEMP, 

KPOSTF, 

KPAT, 

KEIGEN 

KP0ST2 , 

KP0ST3 , 

KP0ST4, 

KP0ST7 , 

KPOST8 , 

KP0ST9 , 
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The common block BOLH contains the logical variable names of program option 
cards . 


COMMON /BOLH/ 

BIDE, 

1 

SUPER, 

2 

POST, 

3 

PLYORD , 

4 

PERIOD, 

5 

RANG, 

6 

SPARS, 

7 

CYLIND, 

8 

EL7 , 

9 

MODLIN 


DYNA, 

PRTOUT , 

WINDML, 

SOLID, 

UNSY, 

PROFIL, 

PROT, 

PREP, 

MODES , 

HOLLOW, 

DISPL, 

DAMP, 

ACCEL, 

TEM, 

VELO, 

PRES, 

PATRAN, 

MHOST, 

EZRD, 

ECHO, 

ANGULA, 

CFORCE, 

PANEL, 

LINEAR, 

RESTAR, 

HEAT, 

MODEL, 

FABRIC, 

BUCK, 

INTER, 

BMODEL, 

EL3 , 

ELIO, 

EL11 , 

EL75 , 

TRANS, 


DYNA 

WINDML 

SOLID 

UNSY 

PROFIL 

MODES 

DISPL 

DAMP 

PLYORD 

ACCEL 

TEM 

VELOL 

PRES 


Flag for dynamic analysis using direct time integration. 
Flag for the HPLATE model option. 

Flag for the HSOLID model option. 

Flag for an unsymmetrical ply lay-up. 

Flag for the PROFILE program option card. 

Flag for modal analysis. 

Flag for inputing initial displacements. 

Flag for damping. 

Flag for symmetrical ply order. 

Flag for inputing initial accelerations. 

Flag for temperature loading. 

Flag for inputing initial velocities. 

Flag for pressure loading. 
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PERIOD 

Flag 

for 

ECHO 

Flag 

for 

ANGOLA 

Flag 

for 

C FORCE 

Flag 

for 

PANEL 

Flag 

for 

SPARS 

Flag 

for 

RESTAR 

Flag 

for 

MODEL 

Flag 

for 

FABRIC 

Flag 

for 

BUCK 

Flag 

for 

INTER 

Flag 

for 

BMODEL 

Flag 

for 

EL3 

Flag 

for 

EL7 

Flag 

for 

ELll 

Flag 

for 

EL75 

Flag 

for 

TRANS 

Flag 

for 

MODLIN 

Flag 

for 
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periodic loading, 
echoing input file, 
centrifugal loading, 
nodal loads. 

initiating the PANEL option. 

generating spar elements. 

initiating a restart run. 

creating a PATRAN neutral file. 

initiating an analysis with fabrication. 

buckling analysis. 

setting an interface. 

the H3DS0LID model option. 

plane stress element. 

3-D solid element, 
plane strain element, 
plate element. 

inputing nodal transformations, 
reading in finite element model. 


The common block INTPH contains the reference temperature. 


COMMON /INTPH/ TREF, BFT, BFC , BCS , FDC, 

1 FDK 
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TREF 

The 

COMMON 

1 

2 

3 

4 

5 

6 
7 

NMECHC 

NTHERC 

LINC 

MSTART 

MITER 

NTISTP 

NCFOR 

TOL 

NEIGV 

NSUBD 

INCREG 

NBC 

MHITER 

Chapter 


Reference temperature 


common block LOADD contains the load data. 


/LOADD/ NMECHC , 

NTHERC , 

LINC, 

MSTART, 

MITER, 

NTISTP, 

NCFOR, 

TOL, 

NEIGV, 

NSUBD, 

INCREG, 

NBC, 

MHITER, 

RESID, 

AESID, 

RERMSQ, 

AERMSQ, 

NPRES , 

NTEMP, 

NTR, 

NBCON , 

NTCON, 

NUMPLY, 

NOPERD, 

DAMPMS, 

DAMPST , 

NDIS, 

NVEL, 

NACC, 

INCDYN 




Number of mechanical cycles. 

Number of thermal cycles . 

Number of load increments between load steps. 

Number of increments to be accumulated before a restart file is 
written. 

Maximum number of allowable iterations for global convergence. 
Number of load steps. 

Number of nodal loads. 

Global tolerence. 

Number of eigenvalues to be extracted. 


Number of subspace dimensions. 

Number of increments between eigenvalue extraction. 
Number of boundary conditions. 

Allowable number of iterations in MHOST. 


4 
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RESID 

NPRES 

NTEMP 

NTR 

NBCON 

NTCON 

NUMPLY 

NOPER 

DAMPMS 

DAMPS! 

NDIS 

NVEL 

NACC 

INCDYN 
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Allowable residual error in MHOST. 

Number of edge loads . 

Number of temperature data sets . 

Number of transformation data sets. 

Scratch variable for boundary conditions. 

Scratch variable for nodal transformations . 

Number of plies. 

Number of periodic loading data sets. 

Damping coefficient applied to the mass matrix. 

Damping coefficient applied to the stiffness matrix. 

Number of initial displacement data sets. 

Number of initial velocity data sets. 

Number of initial acceleration data sets. 

Number of accumulated load increments before material updating. 


The common block ETABLE contains the element data. 


COMMON /ETABLE/ NELCRD , NELNFR, NELNOD, NELSTR , 

1 NELCHR , NELCMP , NDI 


NELCRD Number of coordinates per node. 

NELNFR Number of degrees -of -freedom per node. 

NELNOD Number of nodes per element. 

NELSTR Number of generalized stresses. 
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NELCHR Number of element material properties 
NELCMP Number of stress and/or strain components per node. 
NDI Number of direct stress components per node. 


The common blocks HINDS . SOLDYN . and VALUEH contain model data. Note that 
the common block HINDS is used in the HPLATE model option and the common block 
SOLDYN is used in the HSOLID model option. The common block VALTJEH is used in 
both. 


COMMON 

/HINDS/ 

ALPHA(3 , 3) , 

NGRID , 

NELT, 

NPAST , 

1 


THMAX, 

NXSECT , 

NS PAR, 

NXSPAR, 

2 


NYSPC , 

NYPTS , 

NNQ, 

NNX, 

3 


LSECT(2) , 

NPTOT , 

NETOT , 

HOLD (10) 

4 


MAXPLY, 

MNNX, 

MNNSECT , 

MMWSCT, 

5 


MNXSPC1 , 

MNXSPC , 

MKNT, 

NKNT, 

6 


NLEPM , 

MNSPDES , 

MK, 

MMSCAN, 

7 


SPAMTH, 

PLYMTH , 

MNOPLY , 

NPTOT1 , 

8 

COMMON 

/SOLDYN/ 

MAXNP 

SPMIN, 

SPMAX, 

NGRD, 

MNOPLY, 

1 

LKDES , 

LMAX, 

MAXPLY, 

PLYMTH 

COMMON 

/VALUEH/ 

NP, 

NE, 

NDES , 

IMAT, 

1 


KSMF , 

iu, 

JU, 

XBEG , 

2 


XEND, 

NMODE(3) , 

NOPLY, 

NSECT, 

3 


MESH, 

NGROUP , 

NCOREG , 

IGRD 


NELT 

Scratch variable. 


NPAST 

Scratch variable. 


NXSECT 

Number of input cross sections. 

NS PAR 

Number of spars . 


NXSPAR 

Number of cross sections 

containing spars. 

NYSPC 

Number of elements along 

the Y-axis. 
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NYPTS Number of nodes along the Y-axis. 

NNQ NYSPC + 2 
NNX NXSCET + 1 

LSECT Number of input sections on the bottom and top surfaces . 

NPTOT Number of nodes in the shell portion of a f. e. model generated in the 
HPLATE model option. 

NETOT Total number of elements in f. e. model. 

MAXPLY Maximum number of plies. 

MNNX The maximum of NNX and NXSPAR. 

MNNSECT The maximum of LSCET(l) and LSECT(2). 

MMWSCT The maximum value in the array MWSCT . 

MNXSPC1 MNXSPC + 1 

MNXSPC Maximum value in the array NXSPC plus 1. 

MNSPDES Maximum value in the array NSPDES. 

SPAMTH Minimum spar thickness. 

PLYMTH Minimum ply thickness. 

MNOPLY Maximum number of plies required at any one node. 

MAXNP Total number of nodes . 

SPMIN Minimum ply thickness. 

SPMAX Maximum wall thickness. 

LMAX Number of plies specified for one -half of an unsymmetrical ply lay-up. 

NP Number of nodes . 

NE Number of elements. 

NDES Number of ply descriptions. 
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IU 

JU 

XBEG 

XEND 

NSECT 

IGRD 

The 

COMMON 

1 

ID 

IOPT 

NPRT 

NPRTS 

NPRTP 

NPPLY 

TIMEPN 

TIMEMS 

TIMEMB 


Chapter 


Number of nodes along X-axis. 

Number of nodes along Y-axis. 

Initial x coordinate of structure. 

Final x coordinate of structure. 

Number of cross sections. 

Grid point input options. 

common block PRINH contains the output control information. 

/PRINH/ ID(37 , 5) , I0PT(48,30) , NPRT(2,10) , 

NPRTS (2 , 10) , NPRTP (2 , 10) , NPPLY (2 , 10) , 

TIMEPN (10) , TIMEMS (10) , TIMEMB (10) 

Array containing the header for the 5 title lines. 

Array storing the program options. 

Displacement output sets. 

Ply stress output sets. 

Ply property output sets. 

Ply property output sets. 

Array containing times at which PATRAN result files are to written. 

Array containing times at which eigenvalues are to be extracted. 

Array containing times at which a critical buckling load is to be 
extracted. 
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The common block LCOUNT contains the following variable. 

COMMON /LCOUNT/ LCN1, LCN2, LCN3, LCN4, 

1 LCN5 , LCN6 , LCN7, LCN8 

LCNl Beginning of the temporary storage. 


Chapter 4 


March, 1992 


STOR23 



HITCAN Programmer's Manual - Version 1.0 

CHAPTER 5 

FILE SYSTEM 

The following is a summary of the FORTRAN I/O units used in HITCAN. The 
FORTRAN units are assigned in the subroutine HINIT and reside in the common 
blocks HUNITS and CMETCN . The unit numbers of these files and the variables 
they contain are described below. 

Note that additional files, not listed here, are used in METCAN, COBSTRAN, 
and MHOST. These files are assigned in the subroutines BLOCK, CINIT, and 
BLOCKM. The subroutine BLOCK resides in the METCAN library module, CINIT is 
in the COBSTRAN library module, and the BLOCKM resides in the MHOST library 
module. These files are integral to the library modules in which they are 
used and are not described here. For further information on these files see 
References 4 and 5 . 


FORTRAN 


Code 

Unit No . 

File Contents 

KREAD 

5 

Input (read) unit 

KOUT 

6 

Output (write) unit 

KPAT 

18 

PATRAN neutral file containing the finite element model 

KPOSTF 

19 

Postprocessing file from MHOST 

KSTR 

21 

Contains fiber, matrix, and interface microstresses 
(FSIGNN, MSIGNN, and DSIGNN) and micros tr ess rates 
( FDOTN , MDOTN, and DDOTN) 

KLAM 

22 

Contains flags for material failure LAGF, FLAGM, and FLAGD) 
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FORTRAN 


Code 

Unit No . 

File Contents 

KSIG 

24 

Contains cumulative ply stresses and ply strains (SIGLC and 
EPSPLC) 

KMHOST 

27 

Input file for MHOST from HITCAN 

KP0ST2 

32 

METCAN output file. Contains the current time and fiber 
material properties (TIME and PFN) 

KP0ST3 

33 

METCAN output file. Contains the current time and matrix 
material properties (TIME and PMN) 

KPOST4 

34 

METCAN output file. Contains the current time and interface 
material properties (TIME and PDN) 

KP0ST5 

35 

METCAN output file. Contains the current time and ply 
properties (TIME and PPN) 

KP0ST6 

36 

METCAN output file. Contains the incremental fiber 
microstresses (FS1GN) 

KPOST7 

37 

METCAN output file. Contains the incremental matrix 
microstresses (MSIGN) 

KP0ST8 

38 

METCAN output file. Contains the incremental interface 
microstresses (DSIGN) 

KPOST9 

39 

METCAN output file. Contains the time, incremental ply 
stresses, incremental ply strains, incremental composite 
stresses, and incremental composite strains (TIME, STRNS1, 
STRNS2, STRNS3 , LSIGN, EPSPL, and CSTRES) 

KPOSTO 

40 

METCAN output file. Contains the current composite 
properties 

KBULK 

48 

Scratch file for HITCAN 

KREFL 

49 

For input dataset reflection 

KCBST 

51 

Input file for COBSTRAN 

KCBFE 

54 

Post file from COBSTRAN containing the finite element model 
of a solid structure, using the eight node solid element 

KMETCN 

Chapter 

67 

5 

Input file for METCAN 
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Code 

KSTART 

KDBANK 

KCBFE 

KNODAL 

KTRAN 

KDISP 

KTRANF 

KMHFOR 

KCBOUT 


FORTRAN 

Unit No . File Contents 

68 HITCAN output file containing restart Information 

70 Databank unit. 

75 Post file from COBSTRAN containing the ply lay-up at each 
node. Used for a solid structure modeled with the eight 
node solid element 

76 PATRAN results file containing the nodal displacements 

77 PATRAN results file containing the nodal stresses 

78 MHOST output file containing the incremental displacements 
(DISP) and reactions (REACF) 

89 MHOST output file containing the global to local 
transformation matrix 

93 MHOST output file containing the incremental stress 

resultants (NSB, MSB, and QSB) 

98 Output file from COBSTRAN 
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CHAPTER 6 

HITCAN DICTIONARY 


ACCEL 
ACCELI 
ACELIN 
AES ID 
AERMSQ 

ALPH (1,1) 

ALPH (2,1) 

ALPH (3,1) 

ALPHA 

ANG 

ANGULA 

ANGVEL 

AS 

ANVDIF 

BS 

BMODEL 

BUCK 

Chapter 6 


LOGI Flag for inputing initial accelerations 

REAL Array containing the initial nodal accelerations 

REAL Array containing the initial nodal accelerations 

REAL Maximum allowable absolute error in residuals 

REAL Maximum allowable relative error in the root mean square 

of the strain energy 

REAL Coefficient of thermal expansion in the x- direction at 
node I 

REAL Coefficient of thermal expansion in the y- direction at 
node I 

REAL Coefficient of thermal expansion in the z -direction at 
node I 

REAL Matrix of direction angles between local coordinates and 
global coordinates 

REAL Ply orientation angle 

LOGI Flag for centrifugal loading 

REAL Array containing the angular velocity at each time step 

REAL Axial stiffness matrix 

REAL Incremental angular velocity 

REAL Bending stiffness matrix 

LOGI Flag for the H3DS0LID model option 

LOGI Flag for buckling analysis 
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CFORCE 
CODES (1) 
CODES (2) 
CODES (3) 
CODES (4) 
CODES (S) 
CODES (6) 
CONV 
COFORC 
CS 

CSTRES 

DAMP 

DAMPMS 

DAMPST 

DDOTO 

DDOTN 

DDOTNN 

DDOTI 

DEN 

DFH 

DISP 

DISPI 

DISPIN 
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LOGI Flag for nodal loads 

CHAR Type of ply fiber 

CHAR Type of ply matrix 

REAL Ply thickness 

REAL Void volume ratio 

REAL Fiber volume ratio 

REAL Ply orientation angle (degrees) 

LOGI Flag for global convergence 

REAL Array containing the incremental loads 

REAL Coupling stiffness matrix 

REAL Incremental composite stresses 

LOGI Flag for damping 

REAL Coefficient for the mass matrix for Rayleigh damping 

REAL Coefficient for the stiffness matrix for Rayleigh damping 

REAL Array of initial interface microstress rates 

REAL Array of incremental interface microstress rates 

REAL Array of cumulative interface microstress rates 

REAL Array of interface microstress rates read from the 

databank 

REAL Array of nodal densities 

REAL Array of fiber bundle diameters 

REAL Array of nodal displacements 

REAL Array containing the initial nodal displacements 

REAL Array containing the initial nodal displacements 

March, 1992 


DICT2 



DISPL 

LOGI 

DISVH 

REAL 

DPI 

REAL 

DSIGO 

REAL 

DSIGN 

REAL 

DSIGNN 

REAL 

DSIGI 

REAL 

DTI 

REAL 

DTVCI 

REAL 

DYNA 

LOGI 

ECHO 

LOGI 

EL3 

LOGI 

EL7 

LOGI 

EL11 

LOGI 

EL75 

LOGI 

EPSDC 

REAL 

EPSFC 

REAL 

EPSLC 

REAL 

EPSMC 

REAL 

EPSPL 

REAL 

FABRIC 

LOGI 

FPI 

REAL 

FDOTO 

REAL 

FDOTI 

REAL 

Chapter 6 
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Flag for inputing initial displacements 
Array of incremental displacements 

Array of interface properties read from the databank 

Array of initial interface microstresses 

Array of incremental interface microstresses 

Array of cumulative interface microstresses 

Array of interface microstresses read from the databank 

Incremental elapsed cpu time 

Array of interface exponents read from the databank 

Flag for dynamic analysis using direct time integration 

Flag for echoing input file 

Flag for plane stress element 

Flag for 3-D solid element 

Flag for plane strain element 

Flag for plate element 

Cumulative interface strain 

Cumulative fiber strain 

Cumulative ply strains 

Cumulative matrix strain 

Incremental ply strains 

Flag for initiating an analysis with fabrication 

Array of fiber properties read from the databank 

Array of initial fiber microstress rates 

Array of fiber microstress rates read from the 
databank 
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FDOTN 

REAL 

FDOTNN 

REAL 

FLAGD 

REAL 

FLAGF 

REAL 

FLAGM 

REAL 

FSIGO 

REAL 

FSIGI 

REAL 

FSIGN 

REAL 

FSIGNN 

REAL 

FTVCI 

REAL 

GNP 

REAL 

GRIDPl 

REAL 

GRIDP2 

REAL 

ID 

CHAR 

IGRD 

INT 

INCDYN 

INT 

INCR 

INT 

INCREG 

INT 

INTER 

LOG I 

IOPT 

INT 

IRAMP 

INT 

IU 

INT 

Chapter 6 



Array of incremental fiber microstress rates 

Array of cumulative initial fiber microstress 
rates 

Array of flags for interface property failure 

Array of flags for fiber property failure 

Array of flags for matrix property failure 

Array of initial fiber microstresses 

Array of fiber microstresses read from the databank 

Array of incremental fiber microstresses 

Array of cumulative fiber microstresses 

Array of fiber exponents read from the databank 

Array of nodal coordinates and thicknesses 

First point of the vector required to define the axis of 
rotation 

Second point of the vector required to define the axis of 
rotation 

Array of five title lines 
Grid point input options 

Number of increments per load increment in the direct time 
integration 

Current Increment 

Increment number at which eigenvalues are to extracted 
Flag for setting an interface 
Array storing the program options. 

Current ramp 

Number of nodes along x-axis 
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JU 

INT 

Number of nodes along y-axis 

KAD 

IN 

Used to store shell grid point numbers corresponding to 
the spar internal nodes 

KEI 

INT 

Array of element connectivities 

KF 

REAL 

Apparent fiber volume ratio 

KFB 

REAL 

Actual fiber volume ratio 

KM 

REAL 

Apparent matrix volume ratio 

KMB 

REAL 

Actual matrix volume ratio 

KV 

REAL 

Void volume ratio 

L 

INT 

Number of points in interpolation function 

LINEAR 

LOG I 

Flag for a linear analysis 

LCN1 

INT 

Beginning of the temporary storage. 

LDALOW 

INT 

Number of allowable load steps in METCAN 

LDITER 

INT 

Number of iterations performed in METCAN 

LINC 

INT 

Number of load increments 

LMAX 

INT 

Number of plies specified for opposite half of an 
unsymmetrical lay-up 

LSECT(l) 

INT 

Number of input cross sections on the top surface (HPLATE 
model option) 

LSECT(2) 

INT 

Number of input cross sections on the bottom surface 
(HPLATE model option) 

LSECT 

INT 

Number of input sections of each input plane (S3DS0LID 
model option) 

LSIGN 

REAL 

Incremental ply stresses 

LSMAX 

INT 

Maximum value in the array LSECT 

LSP 

INT 

Location of the spar on the top and bottom surfaces from 
the left end 

Chapter 6 


March, 1992 


DICT5 


MAXNP 

MAXPLY 

HDES 

MDOTO 

MDOTI 

MDOTN 

MDOTNN 

MESH 

MHITER 

MITER 

MMAX 

MMWSCT 

MNNX 

MNNSECT 

MNSPDES 

MNXSPC 

MNXSPC1 

MNOPLY 

MODEL 

MODES 

MODLIN 

MSH 

MPLY 
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INT N umb er of nodal points including both shell and spar nodes 

INT Number of plies specified for half- thickness at point of 

maximum thickness 

INT Array containing ply lay-up for each node 

REAL Array of initial matrix microstress rates 

REAL Array of matrix microstress rates read from the databank 

PEA T. Array of incremental initial matrix microstress rates 

REAL Array of cumulative initial matrix microstress rates 

INT Modeling technique. No longer applicable 

INT Maximum number of iterations allowed In MHOST 

INT Maximum number of iterations per load increment 

INT Maximum value in the array MSECT 

INT The maximum value in the array MWSCT 

INT The maximum of NNX and NXSPAR 

INT The maximum of LSCET(l) and LSECT(2) 

INT Maximum value in the array NSPDES 

INT Maximum value in the array NXSPC plus 1 

INT MNXSPC + 1 

INT Maximum number of plies at a node 

LOGI Flag for creating a PATRAN neutral file 

LOGI Flag for modal analysis 

LOGI Flag for reading in finite element model 
REAL Applied moments at nodal points 

INT Ply order using ply designation numbers for one half of 
the thickness starting at the bottom surface 
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MPI 
MSIGO 
MSIGI 
MSIGN 
MSIGNN 
MSECT „ 
MSTART 
MTV Cl 
MWSCT (J , 1) 

MWSCT (J, 2) 

N 

NACC 

NACCEL(l) 

NACCEL(2) 

NACCEL(3) 

NBC 

NBOUND(l) 
NBOUND ( 2 ) 
NBOUND ( 3 ) 
NBOUND (4) 
NCFOR 
NCFDIR 
NCFNOD 
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REAL Array of matrix properties read from the databank 

REAL Array of initial matrix microstresses 

REAL Array of matrix microstresses read from the databank 

REAL Array of incremental matrix microstresses 

REAL Array of cumulative matrix microstresses 

INT Number of input points for each cross section 

INT Write a restart file after this many load increments 

REAL Array of matrix exponents read from the databank 

INT Number of grid points input on the top surface for each 
cross section J 

INT Number of grid points input on the bottom surface for each 
cross section J 

INT Number of points for interpolation 

INT Number of initial acceleration data sets 

INT Beginning node in an initial acceleration data set 

INT Ending node in an initial acceleration data set 

INT Increment 

INT Number of boundary conditions 

INT Beginning node in a boundary data set 
INT Ending node in a boundary data set 

INT Increment 

INT Degree* of -freedom which Is fixed 

INT Number of concentrated nodal forces 

INT Direction of the concentrated nodal force 

INT Array of nodes having a concentrated nodal force 
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NCFVAL 


REAL 


NDES 

INT 

NDI 

INT 

NDIS 

INT 

NDISPL(l) 

INT 

NDISPL(2) 

INT 

NDISPL(3) 

INT 

NE 

INT 

NEIGV 

INT 

NELCHR 

INT 

NELCMP 

INT 

NELCRD 

INT 

NELNFR 

INT 

NELNOD 

INT 

NELSTR 

INT 

NESPAR 

INT 

NETOT 

INT 

NETT 

structure 

INT 

NEXX 

INT 

NEYY 

INT 

NF 

INT 

NIPL 
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Magnitude of the concentrated force at each time step 
Number of ply designations 

Number of direct stress components per node 
Number of initial displacement data sets 
Beginning node in an initial displacement data set 
Ending node in an initial displacement data set 
Increment 

Number of elements 

Number of modes to be extracted 

Number of element material properties 

Number of stress and/or strain components per node 

Number of coordinates per node 

Number of degrees -of- freedom per node 

Number of nodes per element 

Number of generalized stresses per node 

Number of spar elements 

Total number of elements 

Number of elements through the thickness of a solid 
modeled by brick elements 

Number of elements along the X-axis of a solid structure 
modeled by brick elements 

Number of elements along the Y-axis of a solid structure 
modeled by brick elements 

Array containing the number of fibers per bundle for each 
ply 

Number of input planes 
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NLEP 

INT 

Array containing the number of plies at each node 

NLS 

INT 

Number of load steps in METCAN 

NMECHC 

INT 

Number of mechanical cycles 

NNQ 

INT 

NYSPC + 2 

NNX 

INT 

NXSECT + 1 

NOPERD 

INT 

Number of periodic loadings 

NORAMP 

INT 

NTISTP - 1 

NOSC 

INT 

Number of output sections 

NOSC1 

INT 

NOSC + 1 

NPERD(l) 

INT 

Beginning node in a periodic loading data set 

NPERD ( 2 ) 

INT 

Ending node in a periodic loading data set 

NPERD(3) 

INT 

Increment 

NPERD(4) 

INT 

Loading type 

NPERD(5) 

INT 

Degree -of- freedom the load is applied to 

NP 

INT 

Number of nodes 

NPAST 

INT 

Last node of each SHELL/SPAR cycle 

NPLSET 

INT 

Number of data sets describing the ply lay-up 

NPLY 

INT 

Ply order using ply designation numbers for one half of 
the thickness starting at the top surface, except this 
order is for the opposite half of an unsymmetrical lay-up 

NPNOD(l) 

INT 

Beginning node in a ply lay-up data set 

NPNOD ( 2 ) 

INT 

Ending node in a ply lay-up data set 

NPNOD ( 3 ) 

INT 

Increment 

NPNOD (4) 

INT 

Number of plies in the lay-up 

NPRT (1,1) 

INT 

Initial node for displacement output for set 
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NPRT (2,1) 

INT 

NPRTS(l.I) 

INT 

NPRTS (2,1) 

INT 

NPRTP (1,1) 

INT 

NPRTP(2 , I) 

INT 

NPPLY (1,1) 

INT 

NPPLY(2 ,1) 

INT 

NPRES 

INT 

NPRESS(l) 

INT 

NPRESS(2) 

INT 

NPRESS(3) 

INT 

NPPLY(2 ,1) 

INT 

NPTOT 

INT 

NSECT 

INT 

NSH 

REAL 

NSPAR 

INT 

NSPDES 

INT 

NSUBD 

INT 

NTEMP 

INT 

NTHERC 

INT 

NTISTP 

INT 

NTNOD(l) 

INT 

NTNOD ( 2 ) 

INT 

NTN0D(3) 

INT 
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Final node for displacement output for set I 
Initial node for ply stress output for set I 
Final node for ply stress output for set I 
Initial node for ply property output for set 
Final node for ply property output for set I 
Initial ply for ply property output for set I 
Final ply for ply property output for set I 
Number of edge loads 
Beginning element number 
Ending element number 
Surface with pressure loading 
Final ply for ply property output for set I 
Number of shell elements 
Number of cross sections 
Applied membrane loads at nodal points 
Number of spars 
Spar composite designation 

Number of subspace dimensions for eigenvalue extraction 

Number of temperature data sets 

Number of thermal cycles 

Number of time steps 

Beginning element number 

Ending element number 

Increment 
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Number of plies in this data set 
NTR INT Number of coordinate transformation data sets 

NTRANS(l) INT Beginning node in a transformation data set 

NTRANS(2) INT Ending node in a transformation data set 

NTRANS(3) INT Increment 

NTRANS(4) INT Axis about which the nodal coordinates are to be rotated 
NUMPLY INT Number of plies 

NVEL INT Number of initial velocity data sets 

NVELO(l) INT Beginning node in an initial velocity data set 

NVELO ( 2 ) INT Ending node in an initial velocity data set 

NVELO(3) INT Increment 

NWORDT INT Total number of words required In HITCAN for in core 

storage 

NXDIV INT Array containing the number of elements between cross 

sections 

NXSECT INT Number of cross sections 

NXSPAR INT Number of cross sections containing spars 

NXSPC INT Number of elements between cross sections 

NYPTS INT Number of nodal points along Y-axis 

NYSPC INT Number of elements along Y-axis 

PANEL LOGI Flag for initiating the PANEL option 

PC(1) REAL Composite density 

PC (2) Composite compliance matrix 

through REAL C ( 1 , 1 ) ,C(1,2) ,C(1,3) ,C(2,2) ,C(2,3) , 

PC(ll) C(3,3),C(4,4),C(5,5),C(6,6) 

PC(12) REAL Composite thermal expansion coefficient in 11 direction 
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PC (13) 

REAL 

PC (14) 

REAL 

PC(15) 

REAL 

PC(16) 

REAL 

PC(17) 

REAL 

PC (19) 

REAL 

PC(20) 

REAL 

PC(21) 

REAL 

PC(22) 

REAL 

PC(23) 

REAL 

PC (24) 

REAL 

PC(25) 

REAL 

PC(26) 

REAL 

PC (27) 

REAL 

PC(28) 

REAL 

PC(29) 

REAL 

PC(30) 

REAL 

PDIFU 

REAL 

PDIFL 

REAL 

PDN 

REAL 

PERDLD(l) 

REAL 

PERDLD(2) 

REAL 

PERIOD 

LOGI 

PERT(l) 

REAL 
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Composite thermal expansion coefficient in 22 direction 
Composite thermal expansion coefficient in 33 direction 
Composite thermal conductivity in 11 direction 
Composite thermal conductivity in 22 direction 
Composite thermal conductivity in 33 direction 
Composite modulus in 11 direction 
Composite modulus in 22 direction 
Composite modulus in 33 direction 
Composite shear modulus in 23 direction 
Composite shear modulus in 31 direction 
Composite shear modulus in 12 direction 
Composite poisson's ratio in 12 direction 
Composite poisson's ratio in 21 direction 
Composite poisson's ratio in 13 direction 
Composite poisson's ratio in 31 direction 
Composite poisson's ratio in 23 direction 
Composite poisson's ratio in 32 direction 
Incremental pressure loading on the upper surface 
Incremental pressure loading on the lower surface 
Current interface material properties 
Period of the load 
Intensity of the load 
Flag for periodic loading 

Initial thickness (percent of thickness at each grid 
point) 
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PERT(2) 

REAL 

Final thickness (percent of thickness at each grid point) 

PERT(3) 

REAL 

Initial x-coordinate (percent length) 

PERT (4) 

REAL 

Final x-coordinate (percent length) 

PERT(5) 

REAL 

Initial y-coordinate (percent width) 

PERT(6) 

REAL 

Final y-coordinate (percent width) 

PIN(l) 

REAL 

X coordinate 

PIN(2) 

REAL 

Y coordinate 

PIN(3) 

REAL 

ZU (HSOLID model option) 

Z coordinate (HPLATE & S3DS0LID model options) 

PIN (4) 

REAL 

Z coordinate (HPLATE model options) 
Thickness (HPLATE &S3DS0LID model options) 

PIN(5) 

REAL 

ZL 

PIN(6) 

REAL 

Thickness 

PINTER 

REAL 

Thickness of interface as a percent of fiber diameter 

PFN 

REAL 

Current fiber material properties 

PL 

REAL 

Lower nodal pressure 

PLH(l) 

REAL 

Ply void content 

PLH(2) 

REAL 

Ply apparent fiber content 

PLH(3) 

REAL 

Ply actual fiber content 

PLH (4) 

REAL 

Ply apparent matrix content 

PLH(5) 

REAL 

Ply actual matrix content 

PLH(6) 

REAL 

Ply weight density 

PLH(7) 

REAL 

Ply thickness 

PLH(8) 


Not used 

PLH (9) 


Not used 
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PLH(IO) 

REAL 

Distance from bottom of composite to ply centroid 

PLH(ll) 

REAL 

Distance from reference plane to ply centroid 

PLH(12) 

REAL 

Angle form structural axis to composite material axis 

PLH(13) 

REAL 

Angle from ply material axis to composite material axis 

PLH(14) 

REAL 

Angle from ply material axis to composite structural axis 

PLH(15) 

through 

PLH(23) 

REAL 

Ply stress-strain relations 

SC11 , SC12 , SC13 , SC22 , SC23 , SC33 , SC44 , SC55 , SC66 

PLH(24) 

REAL 

Ply coeff. of thermal expansion in 11 direction 

PLH(25) 

REAL 

Ply coeff. of thermal expansion in 22 direction 

PLH(26) 

REAL 

Ply coeff. of thermal expansion in 33 direction 

PLH(27) 

REAL 

Ply heat conductivity in 11 direction 

PLH(28) 

REAL 

Ply heat conductivity in 22 direction 

PLH(29) 

REAL 

Ply heat conductivity in 33 direction 

PLH(30) 

REAL 

Ply heat capacity 

PLH(31) 

REAL 

Ply modulus in 11 direction 

PLH(32) 

REAL 

Ply modulus in 22 direction 

PLH(33) 

REAL 

Ply modulus in 33 direction 

PLH(34) 

REAL 

Ply shear modulus in 23 direction 

PLH(35 ,1) 

REAL 

Ply shear modulus in 13 direction 

PLH(36 , I) 

REAL 

Ply shear modulus in 12 direction 

PLH(37) 

through 

PLH(42) 

REAL 

Ply Poisson's ratios 

NU12 , NU2 1 , NU13 , NU31 , NU23 , NU32 

PLH(43) 

through 

PLH(49) 


Not used 
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PLH(50) 

PLH(51) 

through 

PLH(71) 

PLH(72) 

PLH(73) 

PLH(74) 

PLH(75) 

PLH(76) 

PLH(77) 

PLH(78) 

PLH(79) 

PLH(80) 

PLYMTH 

PLYORD 

PMN 

PPN 

PRES 

PRESUR 

PREVAL 

PROFIL 

PU 

REACF 

RERMSQ 

RES ID 
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REAL Ply temperature 

Not used 

REAL Ply tensile strength in 11 direction 
REAL Ply compressive strength in 11 direction 
REAL Ply tensile strength in 22 direction 
REAL Ply compressive strength in 22 direction 
REAL Ply tensile strength in 33 direction 
REAL Ply compressive strength in 33 direction 

REAL Ply shear strength in 12 direction 

REAL Ply shear strength in 23 direction 

REAL Ply shear strength in 13 direction 

REAL Minimum ply thickness . 

LOGI Flag for symmetrical ply order 

REAL Current matrix material properties 

REAL Current ply material properties from METCAN 

LOGI Flag for pressure loading 

REAL Incremental pressure loading 

REAL Magnitude of the pressure loading as a function of time 
LOGI Flag for the PROFILE program option card 

REAL Upper nodal pressure 

REAL Incremental nodal reactions 

REAL Maximum allowable relative error in the root mean square 
of the displacement 

REAL Maximum allowable relative error in the residuals 
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RESTAR 

LOG I 

Flag for initiating a restart run 


QSH 

REAL 

Applied shear loads at nodal points 


SIGLC 

REAL 

Cumulative ply stresses 


SOLID 

LOGI 

Flag for the HSOLID model option 


SPARS 

LOGI 

Flag for generating spar elements 


SPAMTH 

REAL 

Minimum spar thickness 


SPMIN 

REAL 

Minimum ply thickness 


SPMAX 

REAL 

Maximum wall thickness 


SY(1,I) 

REAL 

Spar y coordinate from left edge 


SY ( 2 , 1 ) 

REAL 

Spar thickness 


TB 

REAL 

Wall thickness 


TEM 

LOGI 

Flag for temperature loading 


TEPR(1 , I) 

REAL 

Temperature on the upper surface 


TEPR(2 , I) 

REAL 

Temperature on the lower surface 


TEPR(3 , I) 

REAL 

Pressure on the upper surface 


TEPR(4, I) 

REAL 

Pressure on either the lower or inner surface 


TEPLY 

REAL 

Array containing the temperature at each at each node 


TIMEE 

REAL 

Current time 


TIMEPN 

REAL 

Array containing times at which PATRAN plot files are 
written 

to 

TIMEMS 

REAL 

Array containing times at which eigenvalues are to be 
extracted 


TIMEMB 

REAL 

Times at which a critical buckling load is to be extracted 

TIMES (1) 

REAL 

Previous time 


TIMES (2) 

REAL 

Current time 
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THERMO 

REAL 

Array containing the temperature of each ply at a given 
node and a given time step 

THRSTN 

REAL 

Array containing the mid-plane strains and the rotations 
due to thermal loading 

THHF 

REAL 

Half- thickness at each node 

TLAY 

REAL 

Incremental nodal temperature 

TNEW 

REAL 

Current ply temperature 

TOL 

REAL 

Maximum allowable global error on the displacements 

TOLD 

REAL 

Ply temperature at the previous load increment 

TOLMET 

REAL 

Allowable tolerence in METCAN on local convergence 

TOTREA 

REAL 

Cumulative nodal reactions 

TPGNP(l.I) 

REAL 

Temperature on the upper surface of node I 

TPGNP(2 , I) 

REAL 

Temperature on the lower surface of node I 

TPGNP (3,1) 

REAL 

Pressure on the upper surface of node I 

TPGNP(4,I) 

REAL 

Pressure on the lower surface of node I 

TRANG 

REAL 

Angle about which the axis is to be rotated 

TRANS 

LOGI 

Flag for inputing nodal transformations 

TT1 

REAL 

Current elapsed cpu time 

TT2 

REAL 

Previous elapsed cpu time 

U 

REAL 

Array of x values for interpolation 

UNSY 

LOGI 

Flag for an unsymmetrical ply lay-up 

V 

REAL 

Array of y values for interpolation 

VELOL 

LOGI 

Flag for inputing initial velocities 

VELOI 

REAL 

Array containing the initial nodal velocity 

VELOIN 

REAL 

Array containing the initial nodal velocity 
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VDOT 

REAL 

HITCAN Programmer's Manual - Version 1.0 
Array of first derivative of interpolation 

WHAT 

CHAR 

Four character name 

WINDML 

LOG I 

Flag for the HPLATE model option 

X 

REAL 

X coordinate of cross section 

X 

REAL 

Array containing x coordinates of cross sections 

XBEG 

REAL 

Initial x coordinate 

XEND 

REAL 

Final x coordinate 

XINT 

REAL 

Dummy array for interpolation 

XOUT 

REAL 

Dummy array for interpolation 

XX 

REAL 

Array for which the order Is reversed 

Y 

REAL 

Array containing y coordinates 

YINT 

REAL 

Dummy array for Interpolation 

YOUT 

REAL 

Dummy array for interpolation 

YY 

REAL 

Array for which order is reversed 

Z 

REAL 

Z mid- thickness coordinate 

ZB 

REAL 

Distance from bottom of composite to ply centroid 

ZGC 

REAL 

Distance from reference plane to ply centroid 

ZL 

REAL 

Z lower coordinate 

ZU 

REAL 

Z upper coordinate 
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APPENDIX A 
EXECUTIVE MODULE 


C 

C 

c 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


HITCAN 


EXECUTIVE 


MODULE 


METCAN: 

HITCAN: 

MHOST: 


A STORAGE ARRAY 

MAXLEN SIZE OF ARRAY A 

B STORAGE ARRAY 

NASIZE SIZE OF ARRAY B 

IWORK STORAGE ARRAY 

ISIZE SIZE OF ARRAY IWORK 


PROGRAM HITCAN 


LOGICAL 


PLYORD , 

ACCEL, 

TEM, 

PRES, 

VELO, 

1 


PROT, 

PREP, 

POST, 

MODES, 

PRTOUT , 

2 


SOLID, 

SUPER, 

WIND ML, 

DYNA, 

SPARS, 

3 


BIDE, 

UNSY, 

PROFIL, 

HOLLOW, 

DISPL, 

4 


DAMP, 

PERIOD, 

PATRAN, 

MHOST, 

EZRD, 

5 


ECHO, 

RANG, 

ANGULA, 

CFORCE , 

PANEL, 

6 


LINEAR, 

RESTAR, 

HEAT, 

MODEL, 

FABRIC, 

7 


CYLIND, 

BUCK, 

INTER, 

BMODEL, 

ELS, 

8 


EL7 , 

ELIO, 

EL11, 

EL75 , 

TRANS, 

9 


MODLIN 





COMMON 

/HUNITS/ 

KDES1 , 

KKEI , 

KBULK, 

KREAD, 

KOUT, 

1 


KCBST, 

KCBMT , 

KPROP , 

KCBFE , 

KDBANK 

2 


KTRAN, 

KPOST, 

KMHFOR, 

KREFL, 

KPLY, 

3 


KPROPI , 

KSTR, 

KLAM, 

KSIG, 

KNODAL 

4 


KDISP, 

KMSTR, 

KTRANF, 

KCBOUT , 

KTEMP , 

5 


KMHOST , 

KSTART, 

KPOSTF , 

KPAT, 

KEIGEN 

COMMON 

/ BOLH / 

BIDE, 

DYNA, 

PRTOUT , 

WINDML, 

SOLID, 

1 


SUPER, 

UNSY, 

PROFIL, 

PROT, 

PREP, 

2 


POST, 

MODES , 

HOLLOW, 

DISPL, 

DAMP, 

3 


PLYORD , 

ACCEL, 

TEM, 

VELO, 

PRES, 

4 


PERIOD, 

PATRAN, 

MHOST, 

EZRD, 

ECHO, 

5 


RANG, 

ANGULA, 

CFORCE , 

PANEL, 

LINEAR 

6 


SPARS, 

RESTAR, 

HEAT, 

MODEL, 

FABRIC 

7 


CYLIND, 

BUCK, 

INTER, 

BMODEL, 

EL3 , 

8 


EL7 , 

ELIO, 

EL11, 

EL75 , 

TRANS, 

9 


MODLIN 





DIMENSION 

B ( 400000 ) 




DIMENSION 

IWORK ( . 

300000 ) 




DIMENSION 

A ( 8000 

) 




DATA 

NASIZE 

/ 400000 

/ 




DATA 

ISIZE 

/ 300000 

/ 




DATA 

MAXLEN 

/ 8000 / 
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C 

CALL HINIT 

CALL REFLECH ( 1 , KREAD , KREFL , KOUT ) 

CALL LOGOH (KOUT) 

CALL OPTNH ( KBULK , KOUT , KREAD ) 

IF (ECHO) CALL REFLECH ( 2 , KREAD , KREFL , KOUT ) 

CALL OUTHIO (KOUT) 

CALL ELTYPE (KOUT, EL3,EL7, ELIO, EL11,EL75) 

IF (SOLID) CALL HSOLID (A,MAXLEN, IWORK, ISIZE,B,NASIZE) 
IF (WINDML) CALL HPLATE ( A,MAXLEN, IWORK, ISIZE, B,NASIZE) 
IF (BMODEL) CALL S 3D SOL (A, MAXLEN, IWORK, ISIZE, B,NASIZE) 
C 

STOP 

END 
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APPENDIX B 

SUBROUTINE HSOLID 


SUBROUTINE HSOLID 

PURPOSE THIS SUBROUTINE CONTROLS THE GENERATION OF A SOLID 

(THRU-THE-THICKNESS) MODEL. IT CALLS THE APPROPRIATE 
SUBROUTINES TO READ IN THE DATA/ TO SETUP THE DYNAMIC 
ALLOCATION/ TO GENERATE THE FINITE ELEMENT MODEL , AND 
PERFORM THE STRUCTURAL ANALYSIS. 


SUBROUTINE HSOLID 
LOGICAL 

1 

2 

3 

4 

5 

6 

7 

8 
9 
+ 

1 

COMMON /HUNITS/ 

1 

2 

3 

4 

5 

6 
7 

COMMON /BOLH/ 

1 

2 

3 

4 

5 

6 

7 

8 
9 
+ 

1 
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(A,MAXLEN, 

IWORK/ISIZE 

PLYORD/ 

ACCEL, 

VELO, 

PROT/ 

MODES, 

PRTOUT / 

WIND ML, 

DYNA, 

UNSY/ 

PROFIL/ 

DAMP, 

PERIOD, 

EZRD/ 

ECHO, 

CFORCE / 

PANEL, 

HEAT, 

MODEL, 

BUCK, 

INTER, 

EL7 / 

ELIO, 

TRANS, 

MODLIN 

KDES1 / 

KKEI, 

KOUT/ 

KCBST, 

KCBFE/ 

KDBANK, 

KMHFOR/ 

KREFL, 

KSTR/ 

KLAM, 

KDISP / 

KMSTR, 

KTEMP , 

KMHOST, 

KPAT/ 

KEIGEN 

BIDE, 

DYNA, 

SOLID/ 

SUPER, 

PROT/ 

PREP, 

HOLLOW/ 

DISPL, 

ACCEL/ 

TEM, 

PERIOD , 

PATRAN, 

ECHO/ 

RANG, 

PANEL/ 

LINEAR, 

HEAT/ 

MODEL, 

BUCK, 

INTER, 

EL7 / 

ELIO, 

TRANS/ 

MODLIN 


NASIZE) 


TEM, 

PRES, 

PREP, 

POST, 

SOLID, 

SUPER, 

SPARS, 

BIDE, 

HOLLOW, 

DISPL, 

PATRAN, 

MHOST, 

RANG, 

ANGULA, 

LINEAR, 

RESTAR , 

FABRIC, 

CYLIND, 

BMODEL, 

EL3, 

EL11, 

EL75 , 

KBULK, 

KREAD, 

KCBMT , 

KPROP , 

KTRAN, 

KPOST , 

KPLY , 

KPROP I , 

KSIG, 

KNODAL, 

KTRANF, 

KCBOUT , 

KSTART , 

KPOSTF, 

PRTOUT, 

WIND ML, 

UNSY, 

PROFIL, 

POST, 

MODES, 

DAMP, 

PLYORD , 

VELO, 

PRES, 

MHOST, 

EZRD, 

ANGULA, 

CFORCE , 

SPARS, 

RESTAR, 

FABRIC, 

CYLIND, 

BMODEL, 

EL3, 

EL11, 

EL75 , 
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COMMON 

1 

2 

3 

4 

5 

6 

7 

8 

COMMON 

1 

2 

3 

COMMON 

1 

2 

3 

4 

5 

6 

7 

8 
9 
+ 

1 

COMMON 

1 

2 

3 

COMMON 

1 

COMMON 

1 

2 

3 

4 

5 

6 
7 


/POINTS/ 

/PINTTS/ 

/IPOINT/ 

/VALUEH/ 

/SOLDYN/ 

/LOADD/ 

/LCOUNT/ 


COMMON 

1 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 


NM1, NM2 , 
NM8, NM9, 
NM15 , NM16, 
NM22, NM23, 
NM29, NM30, 
NM36, NM37, 
NM43, NM44, 
NM50, NM51 , 
NM57 , NM58, 
NMT61, 
NMT65, 
NMT69, 
NMT73, 

ND1, ND2 , 
ND8, ND9, 
ND15, ND16, 
ND22, ND23 , 
ND29, ND30, 
ND36, ND37 , 
ND43 , ND44, 
ND50 , ND51 , 
ND57 , ND58, 
ND64 , ND65, 
ND71, ND72, 
ND78, ND79, 
NP, 

KSMF, 

XEND, 

MESH, 

SPMIN, 

LKDES, 

NMECHC, 

MITER, 

NEIGV, 

MHITER, 

AERMSQ, 

NBCON , 

DAMPMS , 

NACC, 

LCN1, 

LCN5, 

GRIDP1 ( 3 ) , 
IWORK ( 1 ) 
A ( 1 ) 

B ( 1 ) 


NM3 , NM4 , 
NM10, NM11, 
NM17, NM18, 
NM24, NM25, 
NM31, NM32 , 
NM38, NM39, 
NM45, NM46, 
NM52, NM53, 
NM59 , NM60, 
NMT62 , 
NMT66, 
NMT70, 
NMT74, 

ND3, ND4, 
ND10, ND11 , 
ND17 , ND18, 
ND24, ND25, 
ND31, ND32 , 
ND38, ND39, 
ND45 , ND46, 
ND52, ND53 , 
ND59 , ND60, 
ND66, ND67 , 
ND73 , ND74, 
ND80, ND81, 
NE, 

IU, 

NMODE ( 3 ) , 
NGROUP , 
SPMAX, 

LMAX, 

NTHERC, 

NTISTP, 

NSUBD, 

RESID, 

NPRES, 

NTCON , 

DAMPST, 

INCDYN 

LCN2, 

LCN6, 

GRIDP2 (3 ) , 


NM5, NM6, 
NM12 , NM13 , 
NM19, NM20, 
NM26, NM27 , 
NM33, NM34 , 
NM40, NM41 , 
NM47, NM48, 
NM54, NM55 , 
NM61 , NM62 , 
NMT63, 
NMT67, 
NMT71, 
NMT75, 

NDS, ND6, 
ND12, ND13, 
ND19 , ND20, 
ND26, ND27 , 
ND33, ND34, 
ND40, ND41 , 
ND47 , ND48, 
ND54, ND55, 
ND61, ND62, 
ND68, ND69 , 
ND75, ND76 , 
ND82 , ND83 , 
NDES, 

JU, 

NOPLY, 

NCOREG, 

NGRD, 

MAXPLY, 

LINC, 

NCFOR, 

INCREG, 

AES ID, 
NTEMP , 
NUMPLY, 
NDIS, 


NM7 , 

NM14, 

NM21, 

NM28, 

NM35 , 

NM42, 

NM49 , 

NM56, 

NM63 

NMT64 , 

NMT68, 

NMT72 , 

NMT76 

ND7, 

ND14, 

ND21, 

ND28, 

ND35 , 

ND42 , 

ND49 , 

ND56, 

ND63 , 

ND70 , 

ND77 , 

ND84 

IMAT, 

XBEG, 

NSECT, 

IGRD 

MNOPLY , 

PLYMTH 

MSTART, 

TOL, 

NBC, 

RERMSQ, 

NTR, 

NOPERD, 

NVEL, 


LCN3, LCN4, 
LCN7 , LCN8 
HOLD (10) 
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c 


1 

2 


TT1 = 0.0 

CALL GETIME (TT1,TT2,DTI) 

READ IN PARAMETER DATA 

CALL RINPAR ( KREAD , KOUT , NDES , NSECT , IGRD , IU, «JU , MESH, XBEG , XEND , 

IMAT , NXSPAR , NSPAR , NYSPC, ALPHA , LSECT , NIPL , NOSC , NEYY , 
NETT , MAXNP , NETOT , MAXPLY , LMAX , NPLSET , PINTER ) 

SET THE INITIAL POINTERS FOR DYNAMIC ALLOCATION 

CALL SCPEM1 (NSECT) 

READ IN THE FIRST SET OF MODEL DATA 
CALL SREAD1 ( KREAD, NSECT, B (NM1 ) ) 

SET ADDITIONIAL POINTERS FOR DYNAMIC ALLOCATION 

CALL SCPEM2 (UNSY, HOLLOW, PLYORD,TEM, PRES, ANGULA,CFORCE,D I SPL, 


1 

2 


VELO, ACCEL , NGRD , LMAX , MAXPLY , NPM, B (NM1 ) , MMSECT ,N,N1, 
N2 , N3 , N4 , NSTART ) 

READ IN SURFACE GEOMETRY 

CALL SREAD2 (EL3, ELIO, EL11,EL75, KREAD, KOUT, N3, B(NM1) ,B(NM2) , 

1 B ( NM6 ) ) 

READ IN THE DESCRIPTION OF THE PLIES AND THE PLY LAY-UP 

CALL RINPLY ( PLYORD , UNSY , EL3 , ELIO , ELI 1 , EL7 5 , KREAD , KOUT , NDES , LMAX , 
1 MAXPLY , B ( NM4 ) , B ( NM5 ) , B ( NM8 ) , B ( NM9 ) ) 

DETERMINE THE MAXIMUM NO. OF PLIES IN THE MODEL 

CALL DETNPL (PLYORD, UNSY, N3, B(NM1) ,B(NM2) ,B(NM5) ) 


SET POINTERS REQUIRED FOR ANALYSIS 

CALL INITP3 (KOUT, CFORCE, TRANS, MODES, SUPER, DYNA, BUCK, NSTART, NWORD, 
1 NWORDT , NP , NE , MNOPLY ) 

SET POINTERS FOR THE TEMPORARY WORDS 

CALL SCTEMW (KOUT, NP,N,N1,N2,N3,N4, NWORDT) 
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CHECK FOR ENOUGH MEMORY. IF EXCEEDS NASIZE SET IN THE MAIN 
PROGRAM, TERMINATE THE RUN. 

CALL CHCKMS (KOUT, NASIZE, NSTART,NWORD,NWORDT, JSSTOP) 

IF (JSSTOP .EQ. 1) GO TO 10 

READ IN TEMPERATURE AND PRESSURE DATA 

CALL SREAD3 (TEM, PRES, EL3,EL7, ELIO, EL11,EL75,KREAD, KOUT, N3,N4, 

1 B(NM1) ,B(NM28) ,B(NM35) ,B(NM36) ,B(NMT61) ,B(NMT70) , 

2 B(NMT71) ) 

READ IN ROTATIONAL VELOCITY AND CONCENTRATED FORCE DATA 

CALL RINFOR ( ANGULA , CFORCE , KREAD , KOUT , GRIDP 1 , GRIDP2 , B ( NM2 8 ) , 

1 B (NM29 ) , B(NM30 ) , B(NM31 ) , B (NM32 ) ) 

READ IN DATA FOR DYNMIC ANALYSIS 

IF (DYNA) CALL DYNMIN (EL3, EL7, ELIO, ELll,EL75,DISPL, VELO, ACCEL, 

1 PERIOD, KREAD, KOUT, NP, B(NM39) ,B(NM40) , 

2 B (NM41 ) , B (NM42 ) , B(NM43) ) 

READ IN BOUNDARY CONDITIONS 

CALL RIBOUN ( TRANS, KREAD, KOUT, NBC, NTR, NBCON, NTCON, B ( NM34 ) , 

1 B(NM37) ,B(NM38) ) 

READ IN OUTPUT CONTROL DATA 

CALL OUTPTC (MODES, BUCK, KREAD, KOUT) 

CALL GETIME (TT1,TT2,DTI) 

DTI = DTI 
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C 

c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


CALL THE SUBROUTINE PREPSO TO DETERMINE THE SHELL MODEL AND 
READ IN THE INITIAL MICROMECHANICS MATERIAL PROPERTIES 


CALL PREPSO 

1 

2 

3 

4 

5 

6 

7 

8 
9 
+ 


( PLYORD , UNSY , HOLLOW, TEM, PRES , KOUT, KDBANK, KPROPI , KDES1 , 
KGNP , NTISTP , NTEMP , NUMPLY , N,N2,N3,N4, MMSECT , B ( NM1 ) , 

B ( NM2 ) , B ( NM3 ) , B ( NM4 ) , B ( NM5 ) , B ( NM6 ) , B ( NM7 ) , B ( NM8 ) , 
B(NM9) , B(NM10) ,B(NMll) ,B(NM12) ,B(NM13) , B(NM14) , 

B (NM15 ) , B (NM16) , B (NM17 ) , B(NM18) ,B(NM19) ,B(NM20) , 

B ( NM2 1 ) ,B(NM22) ,B(NM23) ,B(NM25) ,B(NM26) ,B(NM28) , 
B(NM55) ,B(NM56) ,B(NM57) ,B(NM58) ,B(NM59) ,B(NM60) , 
B(NM61) ,B(ND43) ,B(ND44) ,B(ND45) ,B(ND46) ,B(ND47) , 
B(ND48) ,B(ND49) ,B(NMT61) ,B(NMT62) ,B (NMT63 ) , B (NMT64 ) , 
B(NMT65) ,B(NMT66) ,B(NMT67) ,B(NMT68) ,B(NMT69) , 

B(NMT70) ,B(NMT71) ,B(NMT72) ,B(NMT73) ) 


WRITE TO OUTPUT FILE THE PLY STACKUP AT EACH NODE 


CALL OTPUTH ( 3 , SOLID ,WINDML, BMODEL, TEM, PRES , EL3 ,EL7 , ELIO , EL11 , 

1 EL7 5 , KBULK , KOUT , MNOPLY , NSECT , NSPAR , NXSPAR , NPTOT , 

2 NETOT , B ( NM2 8 ) ,B(ND43) ,B(ND44) ,B(ND45) ,B(ND46) , 

3 B(ND47) ,B(NM55) ,B(NM56) ,B(ND48),B(ND42) ) 

WRITE TO OUTPUT FILE THE FINITE ELEMENT MODEL 


CALL OTPUTH (7, SOLID, WINDML, BMODEL, TEM, PRES, EL3,EL7, ELIO, EL11, 

1 EL7 5 , KBULK, KOUT , MNOPLY , NSECT , NSPAR, NXSPAR , NPTOT , 

2 NETOT , B ( NM2 8 ) ,B(ND43) ,B(ND44) , B(ND45) ,B(ND46) , 

3 B(ND47) ,B(NM55) , B(NM56) ,B(ND48) ,B(ND42) ) 

CALL THE SUBROUTINE WNEUT TO GENERATE A PATRAN NEUTRAL FILE 


CALL WNEUT (KPAT, 4, NP, NE, B (ND43 ) , B (ND44) ) 
CALL GETIME (TT1 , TT2 , DTI ) 

DT2 = DTI 

IF (MODEL) GO TO 10 
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c 

c 

c 


c 

c 

c 


c 


c 


c 


CALL THE SUBROUTINE NLINER TO PERFORM THE STATIC ANALYSIS 
IF (.NOT. DYNA) THEN 

CALL NLINER (A, MAXLEN, IWORK, ISIZE, NSTART, IMAT,MESH, NP,NP, NE, 

1 MNOPLY , HOLD , GRIDP1 ,GRIDP2 , PINTER, B(NM5) ,B(NM11) , 

2 B(NM12),B(NM13),B(NM14),B(NM15),B(NM16),B(NM17) , 

3 B(NM18),B<NM19),B(NM20),B(NM21),B(NM22) / B(NM23) , 

4 B(NM25) , B(NM28) , B(NM29) ,B(NM30) f B(NM31) , B(NM32 ) , 

5 B(NM34) , B(NM35) , B(NM36) , B(NM37) , B(NM38) , B(ND42 ) , 

6 B(ND43) ,B(ND44) , B(ND45) , B(ND46) , B(ND47) ,B(ND48) , 

7 B(ND49) , B(ND50) , B(ND51) , B(ND52 ) ,B(ND53) ,B(ND54 ) , 

8 B(ND55) ,B(ND56) , B{ND57) , B(ND58) , B(ND59) ,B(ND60) , 

9 B(ND61) ,B(ND62) , B(ND63) ,B(ND64) r B(ND65) ,B(ND66 ) , 

+ B (ND67 ) f B(ND68) , B(ND69 ) , B(ND70) , B (ND71) , B ( ND72 ) , 

1 B(ND73) ,B(ND74) ,B(ND75) / B(ND76) ) 

CALL THE SUBROUTINE DYNMCH TO PERFORM THE DYNAMIC ANALYSIS 


ELSE 

CALL DYNMCH 

1 

2 

3 

4 

5 

6 

7 

8 
9 
+ 

1 

2 
3 

END IF 


( A , MAXLEN , IWORK, ISIZE , NSTART , IMAT , MESH , NP , NP , NE , 
MNOPLY , HOLD , GRIDP1 , GRIDP2 , PINTER, B ( NM5 ) , B ( NM1 1 ) , 
B (NM12 ) , B(NM13) ,B(NM14) ,B(NM15) ,B(NM16) ,B(NM17) , 
B(NM18) ,B(NM19) ,B(NM20) ,B(NM21) ,B(NM22) ,B(NM23) , 
B(NM25) ,B(NM28),B(NM29) ,B(NM30) ,B(NM31) ,B(NM32) , 
B (NM34 ) , B (NM35 ) , B(NM36 ) , B(NM37 ) , B{NM38) , B (NM39 ) , 
B (NM40 ) , B(NM41 ) , B(NM42 ) , B(NM43 ) ,B(ND42 ) ,B(ND43 ) , 
B (ND44 ) , B(ND45 ) ,B (ND46 ) , B(ND47 ) , B (ND48 ) , B (ND49 ) , 
B (ND50 ) , B(ND51 ) , B (ND52 ) , B (ND53 ) , B{NDS4 ) , B(ND55 ) , 
B(ND56) ,B{ND57) ,B(ND58) ,B(ND59) ,B(ND60) ,B(ND61) , 
B (ND62 ) ,B(ND63) ,B(ND64) ,B(ND65) ,B(ND66) ,B(ND67) , 
B(ND68) ,B(ND69),B(ND70) ,B(ND71) ,B(ND72) ,B(ND73) , 
B (ND74 ) , B(ND75 ) , B(ND76) , B(ND77) ,B(ND78) ,B(ND79) , 
B(ND80) ,B(ND81) ) 


CALL GETIME (TT1,TT2 ,DTI) 
WRITE (KOUT, 20 ) DTI, DT2 , TT1 
10 CONTINUE 


20 FORMAT ( / 1 TIME REQUIRED TO : ’ , 

1 'READ IN DATA 

2 'DO PREPROCESSING 

3 'TOTAL TIME FOR THIS RUN 


• ,18X,F10.3, ’ 
' ,18X,F10.3, ' 
' ,18X,F10.3, ' 


SEC. ’ ,/, 20X, 
SEC. ' , /, 20X, 
SEC. ' ) 


RETURN 

END 
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APPENDIX C 

SUBROUTINE HPLATE 


SUBROUTINE HPLATE 

PURPOSE THIS SUBROUTINE CONTROLS THE GENERATION OF BUILT-UP 
STRUCTURES. IT CALLS THE APPROPRIATE SUBROUTINES TO 
READ IN THE DATA, TO SETUP THE DYNAMIC ALLOCATION, 

TO GENERATE THE FINITE ELEMENT MODEL, AND PERFORM THE 
STRUCTURAL ANALYSIS. 


SUBROUTINE HPLATE 

(A,MAXLEN, 

I WORK, ISIZE, 

B , NASIZE ) 


LOGICAL 


PLYORD , 

ACCEL, 

TEM, 

PRES, 

1 


VELO, 

PROT, 

PREP, 

POST, 

2 


MODES, 

PRTOUT, 

SOLID, 

SUPER, 

3 


WIND ML, 

DYNA, 

SPARS, 

BIDE, 

4 


UNSY, 

PROFIL, 

HOLLOW, 

DISPL, 

5 


DAMP, 

PERIOD, 

PATRAN, 

MHOST, 

6 


EZRD, 

ECHO, 

RANG, 

ANGULA, 

7 


CFORCE , 

PANEL, 

LINEAR, 

RESTAR, 

8 


HEAT, 

MODEL, 

FABRIC, 

CYLIND, 

9 


BUCK, 

INTER, 

BMODEL , 

EL3, 

+ 


EL7 , 

ELIO, 

EL11, 

EL75 , 

1 


TRANS, 

MODLIN 



COMMON 

/HUNITS/ 

KDESl , 

KKEI, 

KBULK, 

KREAD, 

1 

KOUT, 

KCBST, 

KCBMT , 

KPROP , 

2 


KCBFE , 

KDBANK, 

KTRAN, 

KPOST , 

3 


KMHFOR, 

KREFL, 

KPLY, 

KPROP I , 

4 


KSTR, 

KLAM, 

KSIG, 

KNODAL, 

5 


KDISP , 

KMSTR, 

KTRANF, 

KCBOUT , 

6 


KTEMP , 

KMHOST, 

KSTART , 

KPOSTF, 

7 


KPAT, 

KEIGEN 



COMMON 

/BOLH/ 

BIDE, 

DYNA, 

PRTOUT, 

WINDML, 

1 

SOLID, 

SUPER, 

UNSY, 

PROFIL, 

2 


PROT, 

PREP, 

POST, 

MODES, 

3 


HOLLOW, 

DISPL, 

DAMP, 

PLYORD, 

4 


ACCEL, 

TEM, 

VELO, 

PRES, 

5 


PERIOD, 

PATRAN, 

MHOST, 

EZRD, 

6 


ECHO, 

RANG, 

ANGULA, 

CFORCE , 

7 


PANEL, 

LINEAR, 

SPARS, 

RESTAR, 

8 


HEAT, 

MODEL, 

FABRIC, 

CYLIND, 

9 


BUCK, 

INTER, 

BMODEL, 

EL3 , 

+ 


EL7 , 

ELIO, 

EL11, 

EL75 , 

1 


TRANS, 

MODLIN 
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C 


C 

C 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


COMMON /LCOUNT/ 

LCN1 , 

LCN2, 

LCN3, 

LCN4, 

LCN5 , 

LCN6, 

LCN7 , 

LCN8 

COMMON /LOADD/ 

NMECHC, 

NTHERC, 

LINC, 

MSTART 

MITER, 

NTISTP, 

NCFOR, 

TOL, 


NEIGV, 

NSUBD, 

INCREG, 

NBC, 

1 

MHITER, 

RESID, 

AESID, 

RERMSQ 


AERMSQ, 

NPRES, 

NTEMP , 

NTR, 

t 

NBCON , 

NTCON , 

NUMPLY, 

NOPERD 

i 

DAMPMS, 

DAMPST, 

NDIS, 

NVEL, 

r 

NACC, 

INCDYN 



DIMENSION 

NPLY(l) , 

GRIDP1 (3) , 

GRIDP2 ( 3 ) 


DIMENSION 

IWORK ( 1 

) 



DIMENSION 

A ( 1 ) 




DIMENSION 

B ( 1 ) 





TT1 =0.0 
MESH = 4 

CALL GETIME (TT1 , TT2 , DTI ) 

READ IN PARAMETER DATA 

CALL RINPAR (KREAD , KOUT, NDES , NXSECT, IGRD, IU, JU,MESH, XBEG, XEND, 

1 IMAT , NXSPAR, NSPAR , NYSPC , ALPHA, LSECT , NIPL , NOSC , NEYY , 

2 NETT , MAXNP , NETOT , MAXPLY , LMAX , NPLSET , PINTER ) 

SET THE INITIAL POINTERS FOR DYNAMIC ALLOCATION 

CALL HCPEM1 ( PLYORD , ANGULA , CFORCE , PRES, PANEL, KOUT, MPLY,NPM) 

READ IN SPAR DESCRIPTIONS 

CALL RINSPA (KREAD, KOUT, NSPAR, NNX, B (NHl ), B (NH2 ), B (NH3 ), B (NH4 ) ) 

READ IN 'FIRST SET OF MODEL DATA 

CALL HREAD1 (KREAD , KOUT, LSECT,MNNSECT, B (NH33 ) ) 

DETERMINE SEVERAL PARAMETERS REQUIRED FOR DYNAMIC ALLOCATION 

CALL DETPAR ( PLYORD , PANEL, KOUT, B (NHl ), B (NH2 ), B (NH4 ), B (NH5 ), B (NH6 ) , 
1 B ( NH2 2 ) , B ( NH3 3 ) ) 
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SET ADDITIONAL POINTERS FOR DYNAMIC ALLbCATION 
CALL HCPEM2 (CFORCE, PRES, TRANS, NSTART) 

READ IN SURFACE GEOMETRY 

CALL HREAD2 (KREAD,KOUT,B(NH33) ,B(NH34) ) 

READ IN THE DESCRIPTION OF THE PLIES AND THE PLY LAY-UP 

CALL RINPLY ( PLYORD , UNSY , EL3 , ELIO , ELI 1 , EL7 5 , KREAD , KOUT , NDES , LMAX , 

1 MAXPLY,B(NH5),B(NH6),B(NH7),NPLY) 

DETERMINE THE MAXIMUM NUMBER OF PLIES AT A NODE 

CALL HDNPLY (B(NH4) ,B(NH6) ,B(NH33) ,B(NH34) ) 

SET POINTERS REQUIRED FOR THE ANALYSIS 

CALL INITP3 (KOUT, CFORCE, TRANS, MODES, SUPER, DYNA, BUCK, NSTART, NWORD, 
1 NWORDT , MAXNP , NETOT , MNOPLY ) 

SET POINTERS FOR THE TEMPORARY WORDS 

CALL HCTEMW ( CFORCE , TRANS , NWORDT ) 

CHECK FOR ENOUGH MEMORY. IF EXCEEDS NASIZE SET IN THE MAIN 
PROGRAM, TERMINATE THE RUN. 

CALL CHCKMS ( KOUT, NAS I ZE, NSTART, NWORD, NWORDT, JSSTOP) 

IF (JSSTOP .EQ. I) GO TO 5 

READ IN LOAD DATA AND BOUNDARY CONDITIONS 

CALL HREAD3 (PROFIL, PLYORD, TEM, PRES, KREAD, KOUT, B(NH8) ,B(NH27) , 

1 B(NH33) ,B(NH36) ,B(NH37) ,B(NH51T) ,B(NH52T) ) 

READ CONCENTRATED FORCES & ROTATIONAL SPEED 

CALL RINFOR (ANGULA, CFORCE, KREAD, KOUT, GRIDP1, GRIDP2, B(NH8) , 

1 B(NH28) ,B(NH29) ,B(NH30) ,B(NH31) ) 
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C 

C 

C 


C 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


READ IN DATA FOR DYNMIC ANALYSIS 

IF (DYNA) CALL DYNMIN (EL3,EL7, ELIO, ELll,EL75,DISPL,VELO, ACCEL, 

! PERIOD, KREAD, KOUT, MAXNP, B ( NH42 ) ,B(NH43) , 

2 B(NH44) ,B(NH45),B(NH46) ) 

READ IN BOUNDARY CONDITIONS 

CALL RIBOUN ( TRANS, KREAD, KOUT, NBC, NTR, NBCON, NTCON, B ( NH9 ) , 

1 B(NH41) , B(NH42 ) ) 

READ IN OUTPUT CONTROL DATA 

CALL OUTPTC (MODES, BUCK, KREAD, KOUT) 

CALL GETIME (TT1,TT2,DTI) 

DTI * DTI 

CALL THE SUBROUTINE PREPSH TO DETERMINE THE SHELL MODEL 

CALL PREPSH ( SPARS, PLYORD, PRES, TEM,PROFIL, PANEL, CYLIND, B ( NH1 ) , 

1 B ( NH2 ) , B ( NH3 ) , B { NH4 ) , B ( NH5 ) , B ( NH6 ) ,B(NH7) ,B(NH8) , 

2 B(NH9) ,B(NH10) ,B(NH11) ,B(NH12) ,B(NH13) ,B(NH14) , 

3 B(NH15) ,B(NH16) ,B(NH17) ,B(NH18) ,B(NH19) ,B(NH20) , 

4 B(NH21) , B ( NH2 2 ) , B ( NH2 3 ) , B ( NH2 5 ) ,B(NH26) ,B(NH33) , 

5 B(NH34) ,B(NH35) ,B(NH36) ,B(ND42) ,B(ND43) ,B(ND44) , 

6 B (ND45 ) , B (ND46 ) , B (ND47 ) , B(ND48 ) ,B(ND49 ) , B (NH36T) , 

7 B(NH37T) ,B(NH38T) ,B(NH39T) ,B(NH40T) ,B(NH41T) , 

8 B(NH42T) ,B(NH43T) ,B(NH44T) ,B(NH45T) ,B(NH46T) , 

9 B(NH47T) ,B(NH48T) ,B(NH49T) ,B(NH50T) ,B(NH51T) , 

+ B(NH52T) ,B(NH53T) ,B(NH54T) ) 

CALL THE SUBROUTINE PREPSP TO DETERMINE THE SPAR MODEL 

IF (NSPAR .NE. 0) THEN 

CALL PREPSP (SPARS, PLYORD, PRES, TEM,PROT, PANEL, KOUT, B(NH1) ,B(NH2) , 

1 B (NH3 ) , B (NH4) , B (NH5 ) , B(NH6) , B(NH7) ,B (NH8 ) , B(NH22 ) , 

2 B (NH23 ) ,B(NH25) ,B(NH26) ,B(ND42),B(ND43) ,B(ND44) , 

3 B (ND45 ) , B(ND46) ,B(ND47) ,B(ND48) ,B(ND49) ,B(NH51T) , 

4 B(NH52T) ,B(NH53T) ,B(NH54T) ) 

END IF 

CALL THE SUBROUTINE WNEUT TO GENERATE A PATRAN NEUTRAL FILE 
CALL WNEUT (KPAT,4,MAXNP,NETOT,B(ND43) ,B(ND44) ) 
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C 

C 

c 


c 

c 

c 


CALL GETIME (TT1 , TT2 ,DTI ) 

DT2 = DTI 
IF (MODEL) GO TO 5 

CALL THE SUBROUTINE NLINER TO PERFORM THE ANALYSIS 

IF (.NOT. DYNA) THEN . 

CALL NLINER ( A , MAXLEN , IWORK , IS I ZE , NSTART , IMAT , MESH , MAXNP , NPTOT , 
NETOT , MNOPLY , HOLD , GRIDP1 , GRIDP2 , PINTER, B ( NH6 ) , 
B(NH10) , B(NH11) ,B(NH12) ,B(NH13) ,B(NH14) ,B(NH15) , 
B(NH16) ,B(NH17) , B(NH18) ,B(NH19) ,B(NH20) , B(NH21) , 

B(NH23) ,B(NH25) ,B(NH8) ,B(NH28),B(NH29) ,B(NH30) , 

B (NH31) ,B(NH9) ,B(NH27) ,B(NH37) ,B(NH41) ,B(NH42) , 
B(ND42) ,B(ND43) ,B(ND44) ,B(ND45) ,B(ND46) ,B(ND47) , 

B (ND48 ) , B(ND49) , B(ND50) ,B(ND51),B(ND52) ,B(ND53) , 

B (ND54 ) , B (ND55 ) , B(ND56) , B (ND57 ) ,B(ND58) , B (ND59 ) , 
B(ND60 ) , B (ND61 ) , B(ND62 ) , B (ND63 ) , B(ND64) , B(ND65 ) , 
B(ND66) ,B(ND67) ,B(ND68) ,B(ND69) ,B(ND70) ,B(ND71) , 

B (ND72 ) , B (ND73 ) ,B(ND74) ,B(ND75) ,B(ND76) ) 

CALL THE SUBROUTINE DYNMCH TO PERFORM THE DYNAMIC ANALYSIS 


ELSE 

CALL DYNMCH 

1 

2 

3 

4 

5 

6 

7 

8 
9 
+ 

1 

2 

3 


( A , MAXLEN , IWORK , I S I ZE , NSTART , IMAT , MESH , MAXNP , NPTOT , 
NETOT , MNOPLY , HOLD , GRIDP1 , GRIDP2 , PINTER, B (NH6) , 
B(NH10) , B(NH11) ,B(NH12) ,B(NH13) ,B(NH14) ,B(NH15) , 
B(NH16 ) , B (NH17 ) , B (NH18 ) , B (NH19 ) , B (NH20) , B (NH21 ) , 
B(NH23) ,B(NH25) ,B(NH8) ,B(NH28) ,B(NH29) ,B(NH30) , 

B ( NH3 1 ) , B ( NH9 ) , B ( NH2 7 ) ,B(NH37) ,B(NH40) ,B(NH41) , 

B (NH42 ) , B(NH43 ) , B (NH44 ) , B (NH45 ) , B (NH46) , B (ND42 ) , 
B(ND43) ,B(ND44) , B (ND45 ) , B (ND46 ) ,B(ND47 ) , B (ND48 ) , 

B (ND49 ) , B (ND50) , B (ND51 ) , B (ND52 ) ,B(ND53) , B (ND54 ) , 

B (ND55 ) , B (ND56 ) ,B(ND57) ,B(ND58) ,B(ND59) ,B(ND60) , 
B(ND61) ,B(ND62) ,B(ND63) ,B(ND64) ,B(ND65) ,B(ND66) , 

B (ND67 ) ,B (ND68 ) , B(ND69 ) , B (ND70) , B(ND71 ) , B (ND72 ) , 
B(ND73) ,B(ND74) ,B(ND75) ,B(ND76) ,B(ND77) ,B(ND78) , 
B(ND79) ,B(ND80) ,B(ND81) ) 
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END IF 
C 

CALL GETIME (TT1,TT2,DTI) 

WRITE (KOUT, 20 ) DTI, DT2, TT1 
5 CONTINUE 
C 

20 FORMAT ( / 1 TIME REQUIRED TO : ’ , 

1 'READ IN DATA 

2 'DO PREPROCESSING 

3 'TOTAL TIME FOR THIS RUN 

C 

RETURN 

END 


' , 18X, F10. 3, ' SEC. ' » / , 20X, 
• ,18X,F10.3, ' SEC. ' , / , 20X, 
■ ,18X,F10.3, ' SEC.') 
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APPENDIX D 

SUBROUTINE S3DSOL 


SUBROUTINE S3DSOL 

PURPOSE THIS SUBROUTINE CONTROLS THE GENERATION OF STRUCTURES 
USING SOLID ELEMENTS. IT CALLS THE APPROPRIATE SUBR. 
TO READ IN THE DATA, TO SETUP THE DYNAMIC ALLOCATION, 
TO GENERATE THE FINITE ELEMENT MODEL, AND PERFORM THE 
STRUCTURAL ANALYSIS. 


SUBROUTINE S3DSOL 

(A, MAXLEN 

, IWORK, ISIZE, B, NASIZE ) 


LOGICAL 


PLYORD, 

ACCEL, 

TEM, 

PRES, 

VELO, 

1 


PROT, 

PREP, 

POST, 

MODES, 

PRTOUT j 

2 


SOLID, 

SUPER, 

WINDML, 

DYNA, 

SPARS, 

3 


BIDE , 

UNSY, 

PROFIL, 

HOLLOW, 

DISPL, 

4 


DAMP, 

PERIOD, 

PATRAN, 

MHOST, 

EZRD, 

5 


ECHO, 

RANG, 

ANGULA, 

CFORCE, 

PANEL, 

6 


LINEAR, 

RESTAR, 

HEAT, 

MODEL, 

FABRIC, 

7 


CYLIND, 

BUCK, 

INTER, 

BMODEL, 

EL3, 

8 


EL7 , 

ELIO, 

EL11, 

EL75, 

TRANS, 

9 


MODLIN 





COMMON 

/HUNITS/ 

KDES1 , 

KKEI , 

KBULK, 

KREAD, 

KOUT, 

1 

KCBST , 

KCBMT , 

KPROP , 

KCBFE , 

KDBANK, 

2 


KTRAN, 

KPOST, 

KMHFOR, 

KREFL, 

KPLY , 

3 


KPROPI , 

KSTR, 

KLAM, 

KSIG, 

KNODAL 

4 


KDISP, 

KMSTR, 

KTRANF, 

KCBOUT , 

KTEMP , 

5 


KMHOST, 

KSTART , 

KPOSTF, 

KPAT, 

KEIGEN 

COMMON 

/BOLH/ 

BIDE, 

DYNA, 

PRTOUT , 

WINDML, 

SOLID, 

1 

SUPER, 

UNSY, 

PROFIL, 

PROT, 

PREP, 

2 


POST, 

MODES, 

HOLLOW, 

DISPL, 

DAMP, 

3 


PLYORD, 

ACCEL, 

TEM, 

VELO, 

PRES, 

4 


PERIOD , 

PATRAN, 

MHOST, 

EZRD, 

ECHO, 

5 


RANG, 

ANGULA, 

CFORCE , 

PANEL, 

LINEAR 

6 


SPARS, 

RE STAR, 

HEAT, 

MODEL, 

FABRIC 

7 


CYLIND, 

BUCK, 

INTER, 

BMODEL, 

EL3, 

8 


EL7 , 

ELIO, 

EL11, 

EL75 , 

TRANS, 

9 


MODLIN 





COMMON 

/INTPH/ 

TREF, 

BFT, 

BFC, 

BCS, 

FDC, 

1 


FDK 





COMMON 

/LOADD/ 

NMECHC, 

NTHERC, 

LINC, 

MSTART, 

MITER, 

1 

NTISTP, 

NCFOR, 

TOL, 

NEIGV, 

NSUBD, 

2 


INCREG, 

NBC, 

MHITER, 

RESID, 

AESID, 

3 


RERMSQ, 

AERMSQ, 

NPRES, 

NTEMP , 

NTR, 

4 


NBCON, 

NTCON , 

NUMPLY, 

NOPERD, 

DAMPMS 

5 


DAMPST, 

NDIS , 

NVEL, 

NACC, 

INCDYN 
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C 

c 

c 


c 

c 

c 


COMMON /POINTB/ 

1 

2 

3 

4 

5 

6 

7 

8 
9 
+ 

COMMON /IPOINT/ 

1 

2 

3 

4 

5 

6 

7 

8 
9 
+ 

1 

2 

3 

4 

5 

6 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 


NBl, 

NB2 , 

NB3, 

NB4, 

NB5, 

NB6, 

NB7 , 

NB8, 

NB9, 

NB10, 

NB11, 

NB12 , 

NB13 , 

NB14, 

NB15, 

NB16, 

NB17, 

NB18, 

NB19, 

NB20, 

HB21, 

NB22 , 

NB23 , 

NB24 , 

NB25, 

NB26, 

NB27 , 

NB28, 

NB29 , 

NB30, 

NB31, 

NB32, 

NB33, 

NB34, 

NB35 

NB36, 

NB37, 

NB38, 

NB39, 

NB40, 

NB41 , 

NB42, 

NB43, 

NB44, 

NB45 

NB46, 

NB47, 

NB48, 

NB49, 

NB50 

NB51, 

NB52, 

NB53 , 

NB54 , 

NB55 

ND1, 

ND2, 

ND3, 

ND4 , 

ND5, 

ND6, 

ND7 , 

ND8, 

ND9, 

ND10 

ND11, 

ND12, 

ND13, 

ND14 , 

ND15 

ND16, 

ND17, 

ND18, 

ND19 r 

ND20 

ND21 , 

ND22, 

ND23, 

ND24 , 

ND25 

ND26, 

ND27 , 

ND28, 

ND29 , 

ND30 

ND31, 

ND32, 

ND33, 

ND34, 

ND35 

ND36, 

ND37, 

ND38 , 

ND39, 

ND40 

ND41, 

ND42, 

ND43, 

ND44, 

ND45 

ND46, 

ND47 , 

ND48, 

ND49, 

ND50 

ND51, 

N052 , 

ND53, 

ND54 , 

ND55 

ND56, 

ND57, 

ND58, 

ND59, 

ND60 

ND61, 

ND62, 

ND63 , 

ND64, 

ND65 

ND66, 

ND67 , 

ND68, 

ND69, 

ND70 

ND71 , 

ND72, 

ND73, 

ND74, 

ND75 

ND76, 

ND77 , 

ND78, 

ND79, 

ND80 

ND81, 

ND82, 

ND83, 

ND84 



GRIDPl ( 3 ) , GRIDP2 ( 3 ) 

IWORK ( 1 ) 

A ( 1 ) 

B ( 1 ) 


READ IN PARAMETER DATA 


TT1 = 0,0 

CALL GETIME (TT1 , TT2 , DTI ) 

CALL RINPAR (KREAD,KOUT, NDES,NSECT, IGRD, IU, JU, MESH, XBEG, XEND , 

1 IMAT , NXSPAR, NSPAR , NYSPC , ALPHA, LSECT , NIPL , NOSC , NEYY , 

2 NETT , MAXNP , NETOT , MAXPLY , LMAX , NPLSET , PINTER ) 

SET THE INITIAL POINTERS FOR THE DYNAMIC ALLOCATION 


CALL BCPEM1 ( NOSC, NIPL, NOSC1 ) 
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READ IN THE FIRST SET OF DATA DESCRIBING THE SURFACE GEOMETRY 

IF (.NOT. MODLIN) CALL BREAD 1 (KREAD, KOUT, NIPL, NOSC, NOSC1, B(NB1) , 

1 B(NB2) ,B(NB3) ) 

SET ADDITIONIAL POINTERS FOR DYNAMIC ALLOCATION 

CALL BCPEM2 (MODLIN, NOSC, NIPL, NOSC1, LSMAX, NEXX,NEYY, NETT, NPXX, 

1 NPYY , NPTT , MAXNP , NETOT , B ( NB 1 ) , B ( NB3 ) ) 

READ IN ADDITIONIAL DATA DESCRIBING THE SURFACE GEOMETRY 

IF (.NOT. MODLIN) CALL BREAD 2 (KREAD,KOUT,NIPL,LSMAX,NOSCl,MMAX, 

1 B(NB3) , B(NB4) ) 

SET ADDITIONIAL POINTERS FOR DYNAMIC ALLOCATION 

CALL BCPEM3 ( PLYORD , UNSY , CFORCE , ANGULA, PRES, TRANS, NDES,NOSC,NIPL, 

1 MAXPLY , LMAX , NCFOR , NT I STP , NPRES , NBC , NOSC1 , LSMAX , MMAX , 

2 MAXNP, NETOT, MNOPLY, NTR, NOPERD, NWORDT, NSTART, B(NB3) , 

3 B (NB4 ) ) 

SET POINTERS REQUIRED FOR ANALYSIS 

CALL INITP3 (ROUT, CFORCE, TRANS, MODES, SUPER, DYNA, BUCK, NSTART, NWORD, 
1 NWORDT, MAXNP, NETOT, MNOPLY) 

CHECK FOR ENOUGH MEMORY. IF EXCEEDS NASIZE SET IN THE MAIN 
PROGRAM, TERMINATE THE RUN. 

CALL CHCKMS (KOUT, NASIZE , NSTART , NWORD , NWORDT, JSSTOP) 

IF (JSSTOP .EQ. 1) GO TO 5 

READ IN FINITE ELEMENT MODEL 

IF (MODLIN) CALL READML (KREAD, MAXNP, NETOT, B(ND43) ,B(ND44) ) 

READ IN SURFACE GEOMETRY 

IF (.NOT. MODLIN) CALL BREAD 3 ( KREAD, KOUT, NIPL, LSMAX, MMAX, B ( NB3 ) , 

1 B (NB4 ) , B (NB5 ) ) 
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C 

C 

C 


C 

C 

C 


C 

C 

C 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 

c 

c 

c 


READ IN THE DESCRIPTION OF THE PLIES 

CALL RINPLY ( PLYORD , UNSY , EL3 , ELIO , ELI 1 , EL7 5 , KREAD , KOUT , NDES , LMAX , 

1 MAXPLY,B(NB7) , B(NB8) ,B(NB9) / B(NB10) ) 

READ IN THE PLY LAY-UP 

IF (MODLIN) CALL RINLAY ( KREAD , KOUT , NPLSET , MNOPLY , MAXNP , B ( NB8 ) , 
l B(NB9) ,B(ND46) ,B(ND47) , B(ND48) ) 

READ IN TEMPERATURE AND PRESSURE DATA 

CALL BREAD4 ( CFORCE, PRES , TEM, KREAD, KOUT ,NIPL, MAXNP, B (NB 6) ,B(NB11) , 
1 B (NB32 ) , B (NB33 ) , B (ND45 ) ) 

READ CONCENTRATED FORCES & ROTATIONAL SPEED 

CALL RINFOR (ANGULA, CFORCE, KREAD, KOUT ,GRI DPI, GRIDP2 , B(NB11 ) , 

1 B(NB28) ,B(NB29) ,B(NB30) , B (NB31 ) ) 

READ IN DATA FOR DYNMIC ANALYSIS 

IF {DYNA) CALL DYNMIN (EL3,EL7 f EL10,ELll / EL75,DISPL / VELO,ACCEL / 

1 PERIOD, KREAD f KOUT, MAXNP, B(NB3 6) ,B(NB37) , 

2 B (NB38 ) , B (NB39 ) , B (NB40 ) ) 

READ IN BOUNDARY CONDITIONS AND TRANSFORMATIONS 

CALL RIBOUN ( TRANS, KREAD, KOUT, NBC,NTR, NBCON, NTCON, B ( NB12 ) ,B(NB34) , 
1 B (NB35 ) ) 

READ IN OUTPUT CONTROL DATA 

CALL OUTPTC (MODES , BUCK, KREAD, KOUT) 

READ IN THE DATABANK 

CALL BANKRH (KDBANK, KOUT, NDES , B (NB8 ) ,B(NB26) ,B{NB25 ) , B(NB27 ) , 

1 B(NB13) ,B(NB14) ,B(NB15) ,B(NB16) ,B(NB17) ,B(NB18) , 

2 B (NB19 ) , B ( NB20 ) ,B(NB21) ,B(NB22) ,B{NB23) ,B(NB24) , 

3 TREF , BFT , BFC , BCS , FDC , FDK) 

CALL GETIME (TT1 , TT2 , DTI ) 

DTI = DTI 
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CREATE THE INPUT DECK FOR COBSTRAN GENERATE FE MODEL 

IF (.NOT. MODLIN) THEN 
MTCODE = 0 

CALL PCOBST ( PLYORD , UNSY , KREAD , KOUT , KCBST , MTCODE , NDES , NIPL , NOSC , 

1 NOSC1 , LSMAX , MMAX , NEYY , NETT , MAXPLY , LMAX , B ( NB1 ) , 

2 B (NB2 ) , B (NB3 ) ,B (NB4) , B(NB5 ) , B(NB7 ) , B(NB8) , B(NB9 ) , 

3 B(NB10)) 

GENERATE THE FINITE ELEMENT MODEL 

CALL COBSTRAN 

INTERPRET NASTRAN BULKDATA 

CALL NASINT (KCBFE, KOUT, NEXX, NEYY, NETT, NPXX , NPYY , MAXNP , NETOT , 

1 B(ND43) ,B(ND44) ) 

CREATE THE INPUT DECK FOR COBSTRAN CREATE PLY LAY-UP 

MTCODE m 1 

CALL PCOBST (PLYORD, UNSY, KREAD, KOUT, KCBST, MTCODE, NDES, NIPL, NOSC, 

1 NOSC1, LSMAX, MMAX, NEYY, NETT, MAXPLY, LMAX, B ( NB1 ) , 

2 B (NB2 ) , B (NB3 ) , B (NB4 ) ,B(NB5) ,B(NB7) ,B(NB8) ,B(NB9) , 

3 B(NB10)) 

GENERATE THE PLY LAY-UP 
CALL COBSTRAN 

INTERPRET BULKDATA TO GET PLY LAY-UP 

CALL METINT ( KCBMT , KOUT , NEXX , NEYY , NPXX , NPYY , MAXNP , MNOPLY , B ( NB8 ) , 
1 B (ND44 ) , B (ND46 ) ,B(ND47) ,B(ND48) ) 

END IF 

OUTPUT MODEL DATA 

CALL OMODEL ( TEM, PRES , KOUT, MAXNP, NETOT, MNOPLY, NTISTP, B (NBll ) , 

1 B (ND43 ) , B (ND44) ,B(ND45) ,B(ND47) ,B(ND48) ) 

CALL THE SUBROUTINE WNEUT TO GENERATE A PATRAN NEUTRAL FILE 

CALL WNEUT (KPAT, 8, MAXNP, NETOT, B (ND43 ), B (ND44) ) 
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c 

c 

c 


c 

c 

c 


c 

c 

c 


DETERMINE NODAL TRANSFORMATION MATRIX 

CALL BNORML (KOUT, MAXNP, NETOT, B(ND43) ,B(ND44) , B(ND75) ) 

CALL GETIME (TT1,TT2 ,DTI ) 

DT2 = DTI 
IF (MODEL) GO TO 5 

CALL THE SUBROUTINE NLINER TO DO THE STATIC ANALYSIS 
IF (.NOT. DYNA) THEN 

CALL NLINER (A,MAXLEN, IWORK, ISIZE,NSTART, IMAT , MESH , MAXNP , MAXNP , 
NETOT , MNOPLY , HOLD , GRIDP1 , GRIDP2 , PINTER , B ( NB8 ) , 

B (NB13 ) ,B(NB14) ,B(NB15) ,B(NB16) ,B(NB17) ,B(NB18) , 

B(NB19) ,B(NB20) ,B(NB21) ,B(NB22) ,B(NB23) ,B(NB24) , 
B(NB25) ,B(NB26) ,B(NB11) ,B(NB28) ,B(NB29) ,B(NB30) , 

B (NB31 ) , B (NB12 ) ,B(NB32) ,B(NB33) ,B(NB34) ,B(NB35) , 
B(ND42) , B (ND43 ) ,B(ND44) , B(ND45) ,B(ND46) , B(ND47) , 
B(ND48) ,B(ND49) ,B(ND50) ,B(ND51) ,B(ND52) ,B(ND53) , 
B(ND54) , B (ND55 ) ,B(ND56) ,B(ND57) ,B(ND58) ,B(ND59) , 
B(ND60) ,B(ND61) ,B(ND62) ,B(ND63) ,B(ND64) ,B(ND65 ) , 
B(ND66) , B(ND67 ) , B(ND68) f B(ND69) ,B(ND70) ,B(ND71) , 

B (ND72 ) , B (ND73 ) , B(ND74 ) , B(ND75 ) , B(ND76) ) 

CALL THE SUBROUTINE DYNMCH TO PERFORM THE DYNAMIC ANALYSIS 


ELSE 

CALL DYNMCH 

1 

2 

3 

4 

5 

6 

7 

8 
9 
+ 

1 

2 

3 

END IF 


(A,MAXLEN, IWORK, ISIZE,NSTART, IMAT, MESH, MAXNP, NPTOT, 
NETOT , MNOPLY , HOLD , GRI DP 1 , GRIDP2 , PINTER , B ( NB8 ) , 

B(NB13),B(NB14) ,B(NB15),B(NB16),B(NB17),B(NB18), 
B(NB19) ,B(NB20) ,B(NB21) ,B(NB22) ,B(NB23) ,B(NB24) , 

B(NB25) , B (NB26) ,B(NB11) ,B(NB28) ,B(NB29) ,B(NB30) , 
B(NB31) , B(NB12) ,B(NB32) ,B(NB33) ,B(NB34) ,B(NB35) , 
B(NB36) ,B(NB37) ,B(NB38) ,B(NB39) ,B (NB40) , B (ND42 ) , 
B(ND43) , B (ND44 ) ,B(ND45) ,B(ND46) ,B(ND47) ,B(ND48) , 
B(ND49) ,B(ND50) ,B(ND51) ,B(ND52) ,B(ND53) ,B(ND54) , 

B (ND55 ) , B(ND56) , B (ND57 ) , B (ND58) ,B(ND59 ) , B (ND60 ) , 

B (ND61) , B (ND62 ) , B (ND63 ) , B(ND64) ,B(ND65 ) ,B(ND66) , 
B(ND67 ) , B (ND68) , B (ND69 ) , B(ND70) , B(ND71) , B(ND72 ) , 

B(ND73),B(ND74) ,B(ND75),B(ND76),B(ND77),B(ND78), 

B(ND79 ) , B (ND80) ,B(ND81) ) 
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CALL GETIME (TT1 , TT2 , DTI ) 

WRITE ( KOUT ,20) DTI, DT2 , TT1 
C 

5 CONTINUE 
C 

20 FORMAT (/’ TIME REQUIRED TO : 1 , 

1 ’READ IN DATA 

2 ’DO PREPROCESSING 

3 ’TOTAL TIME FOR THIS RUN 

C 

RETURN 

END 


’ ,18X,F10.3, ’ SEC. 
’ ,18X,F10.3, ’ SEC. 
’ ,18X,F10.3, ’ SEC. 


tit 2 OX, 
, / , 20X, 

) 
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APPENDIX E 

SUBROUTINE NLINER 


SUBROUTINE NLINER 

PURPOSE PERFORMS THE INCREMENTAL STATIC ANALYSIS 


SUBROUTINE NLINER (A,MAXLEN, IWORK, ISIZE,NSTART, IMAT,MESH,MAXNP, 


1 


NPTOT , NETOT , MNOPLY , HOLD , GRIDP1 , GRIDP2 

, PINTER, 

2 


CODEH,FPI, 

MPI , DPI , FSIGI ,MSIGI 

, DSIGI , FDOTI , 

3 


MDOTI , DDOTI , FTVCI , MTVCI , DTVCI 

,NF, DFH, 

TISTPS, 

4 


ANGVEL , NCFNOD , NCFDIR 

, CFVAL , NBOUND , NPRE S S , 

5 


PREVAL , NTRANS , TRANG , 

KAD , KE I , GNP , TPGNP 

, THHF , 

6 


NLEP,MDES, 

THERMO, AS, 

CS , BS , ALPH , DEN , TLAY , PD I FO f 

7 


PDIFI,DISP,DISVH,MSH 

, NSH, QSH, 

REACF , COFORC , PLH , 

8 


TOLD , TEMP , 

TNEW , E I GNVA , EIGNVE , 

TEMPER , PC , S TRMOM , 

9 


PRE SUR , ETRAN , THSTRN ) 




REAL 


MTVCI , 

MDOTI, 

MSIGI, 

MPI, 

NSH, 

1 


MSH, 

NS, 

MS 



LOGICAL 


PLYORD, 

ACCEL, 

TEM, 

PRES, 

VELO, 

1 


PROT, 

PREP, 

POST, 

MODES, 

PRTOUT , 

2 


SOLID, 

SUPER, 

WINDML, 

DYNA, 

SPARS, 

3 


BIDE, 

UNSY, 

PROFIL, 

HOLLOW, 

DISPL, 

4 


DAMP, 

PERIOD , 

PATRAN, 

MHOST, 

EZRD, 

5 


ECHO, 

RANG, 

ANGULA, 

CFORCE , 

PANEL, 

6 


LINEAR, 

RESTAR, 

HEAT, 

MODEL, 

FABRIC, 

7 


CYLIND, 

BUCK, 

INTER, 

BMODEL, 

EL3, 

8 


EL7, 

ELIO, 

EL11, 

EL75 , 

TRANS, 

9 


MODLIN 





LOGICAL 


CONV 





COMMON 

/HUNITS/ 

KDES1, 

KKEI, 

KBULK, 

KREAD, 

KOUT, 

1 


KCBST, 

KCBMT , 

KPROP , 

KCBFE, 

KDBANK, 

2 


KTRAN, 

KPOST, 

KMHFOR, 

KREFL, 

KPLY, 

3 


KPROPI , 

KSTR, 

KLAM, 

KSIG, 

KNODAL, 

4 


KDISP, 

KMSTR, 

KTRANF, 

KCBOUT , 

KTEMP , 

5 


KMHOST, 

KSTART, 

KPOSTF , 

KPAT, 

KEIGEN 

COMMON 

/CMETCN/ 

KPOSTO , 

KPOST1 , 

KPOST2 , 

KPOST3 , 

KPOST4 , 

1 


KPOST5 , 

KPOST6, 

KPOST 7, 

KPOST8 , 

KPOST9 , 

2 


KMETCN 
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COMMON /BOLH/ 

1 

2 

3 

4 

5 

6 

7 

8 
9 

COMMON /INTPH/ 

1 

COMMON /LOADD/ 

1 

2 

3 

4 

5 

6 
7 

COMMON /ETABLE/ 

1 

DIMENSION 

1 

2 

3 

4 

5 

6 
7 

DIMENSION 

1 

2 

3 

4 

5 

6 
7 

DIMENSION 

1 

2 

3 

4 


BIDE/ 

DYNA, 

PRTOUT 

, WINDML, SOLID, 

SUPER/ 

UNSY, 

PROFIL 

, PROT, PREP, 

POST/ 

MODES , 

HOLLOW 

, DISPL, DAMP, 

PLYORD/ 

ACCEL, 

TEM, 

VELO, PRES, 

PERIOD, 

PATRAN, 

MHOST, 

EZRD, ECHO, 

RANG/ 

ANGULA, 

CFORCE 

, PANEL, LINEAR 

SPARS/ 

RESTAR, 

HEAT, 

MODEL, FABRIC 

CYLIND, 

BUCK, 

INTER, 

BMODEL, EL3, 

EL7 , 

ELIO, 

EL11, 

EL75 , TRANS, 

MODLIN 




TREF/ 

BFT, 

BFC, 

BCS, FDC, 

FDK 




NMECHC/ 

NTHERC, 

LINC, 

MSTART, 

MITER/ 

NTISTP, 

NCFOR, 

TOL, 

NEIGV/ 

NSUBD, 

INCREG 

, NBC, 

MHITER/ 

RESID, 

AESID, 

RERMSQ, 

AERMSQ/ 

NPRES, 

NTEMP, 

NTR, 

NBCON / 

NTCON , 

NUMPLY 

, NOPERD, 

DAMPMS/ 

DAMPST, 

NDIS, 

NVEL, 

NACC/ 

INCDYN 



NELCRD , 

NELNFR, 

NELNOD 

i, NELSTR, 

NELCHR/ 

NELCMP, 

NDI 


CODEH { 6 / 1 ) / 

TISTPS ( 1) , 

TIMES ( 2 ) , 

FPI ( 22 / 1 ) / 

MPI (22, 

i). 

DPI (22,1) , 

FSIGI (6,1) , 

MSIGI (14,1), 

DSIGI ( 10, 1) , 

FDOTI (6,1) , 

MDOTI (14,1), 

DDOTI (10,1), 

FTVCI ( 31 , 1 ) 

, MTVCI (31,1), 

DTVCI (31,1), 

NF ( 1 ) , 

DFH(l) , 


COFORC ( 1 ) , 

ANGVEL ( 1 ) , 

NCFNOD ( 1 ) , 

NCFDIR(l) , 

CFVAL(NCFOR 

.1). 


PDIFI(l) 

THHF ( 1 ) , 

NLEP ( 1 ) 

F 

CS (NELSTR, 1) , 

DEN(l) , 

TLAY ( 5 , 

Dr 

PDIFO(l) , 

GRIDPl (3 ) , 

GRIDP2 ( 3 ) , 

NSH (3,1) , 

QSH ( 2 , 1 ) , 

KAD(l) , 


BS (NELSTR, 1) , 

TRANG ( 1 ) , 

NTRANS(4,1) , 

ETRAN (3,3,1) , 

DISP (NELNFR 

/ 1 ) f 


REACF ( NELNFR , 1 ) , 

DISVH (NELNFR, 1 ) , 


MSH (3,1) , 

AS(NELSTR,1) r 


GNP( 12, 1) 

TOLD (MNOPLY 



TEMP (MNOPLY, 1) , 

TNEW (MNOPLY 



NBOUND (4,1) , 

TEMPER ( 1 ) , 



PC (62,1), 

STRMOM ( NELSTR , 1 ) , 


PRESUR(l) , 

NPRESS (3,1) 

1 


PREVAL ( NPRES , 1 ) 
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c 


c 

c 

c 


DIMENSION 

1 

2 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 


KEI (9,1), TPGNP { 7 , MAXNP , 1 ) , 

PLH (80,1) , THERMO (MNOPLY, MAXNP, 1) , 

ALPH (3,1), MDES (MNOPLY, 1) 
EIGNVA(l), EIGNVE(NELNFR, MAXNP, 1) 

THSTRN (6,1) 

IWORK ( 1 ) 

A ( 1 ) 


DERORO = 0.0 
TOLMET =0.05 
LDALOW = 10 
LSALOW = 5 
ICONVE = 0 
KCODE = 0 
ICODE = 1 
MMODE = 0 
MBUCK = 0 

IF (MODES) MMODE = 1 
IF (BUCK) MBUCK = 1 
MODES = .FALSE. 

BUCK = .FALSE. 

CALL GETIME (TT1,TT2,DTI) 


SUBROUTINE INITSR INITIALIZE THE STRESS RESULTANTS 


CALL INITSR (MAXNP, MSH,NSH,QSH,SIGZ,PDIFO,PDIFI) 


C 

C IF THIS IS A RESTART RUN, READ DATA FROM FILE KREAD 

C 

IF (.NOT. RESTAR) GO TO 5000 
JCODE = 2 

CALL RESTART (JCODE, KREAD, KSTART, KLAM, KS IG,KSTR,NMC, NORAMP, IRAMP, 
1 I NCR , ACCT I M , NLEP , D I SP , TOTRE A , MAXNP , NELNFR ) 

CALL GETIME (TT1,TT2,DTI) 

ACCT = ACCTIM + DTI 
WRITE (KOUT, 2018) DTI, ACCT 
ACCT = ACCT - DTI 
GO TO 6000 
5000 CONTINUE 


C 

C SUBROUTINE INITMS INITIALIZE THE MICRO STRESSES, ETC. 

C 

CALL INITMS (MAXNP, KSTR,KSIG, KLAM, NLEP, NELNFR, DISP,TOTREA) 
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DETERMINE THE THERMAL & MECHANICAL LOADS IN THE SUBROUTINE LOADIN 

CALL LOADIN ( PRES , TEM , ANGULA , CFORCE , 0 , MAXNP ,1,0, NCFOR , NPRES , CFVAL , 

1 COFORC , TPGNP , PDIFO , PDIFI , TREF, TLAY , ANGVEL , ANVDIF , 

2 NTISTP, PREVAL, PRESUR) 

TIMES (1) = TISTPS ( 1 ) 

TIMES (2) * TISTPS ( 1 ) 

TIMEE = TISTPS (1) 


C ITER. OVER EACH NODE & RUN THE MAT’L PROP. MODULE HTMET (METCAN) . 

C SUB. PMETCN SETS UP THE INPUT FILE FOR METCAN. SUB. PLY DETER. 

C THE PLY PROP. ARRAY PLH ( 1 , PLH{ 14, ). THE SUB. PLPROP READS 

C KPOST5 & SETS UP THE PLY PROP. PLH(15, ) , . . . ,PLH(80, ) FOR 

C EACH PLY. CALL THE SUB. PROP TO OBTAIN THE NODAL PROPERTIES. 

C 

TOTTr* s 1 

DO 30 IK = 1, MAXNP 
C PRINT * , ’ NODE # ’ , IK 

NL = NLEP(IK) 

HTHICK = THHF(IK) 

DO 35 IL = 1 , NL 
C PRINT * , * PLY # * , IL 

TOLD (IL, IK) = TREF 
TNEW ( IL, IK) = TREF 
TEMP( IL, IK) = 0.0 
IF (TEM) THEN 

IF (EL3 .OR. EL 7 .OR. EL11) THEN 
TOLD( IL, IK) = TPGNP ( 6, IK, 1) 

TNEW( IL, IK) * TPGNP(6, IK, 1) 

TEMP (IL, IK) = TPGNP (6, IK, 1) - TREF 
ELSE 

CJOE TOLD (IL, IK) = TREF 
TOLD (IL, IK) =0.0 

IF (FABRIC) TOLD ( IL, IK) = THERMO (IL, IK, 1) 

TNEW( IL, IK) = THERMO (IL, IK, 1) 

TEMP ( IL, IK) = THERMO (IL, IK, 1) - TREF 
END IF 
END IF 
35 CONTINUE 
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PO = DABS(PDIFO(IK) ) 

PI = DABS (PDIFI (IK) ) 

NLS = 2 

CALL PMETCN (IK, IREC, JREC, INTER, KSTR, KOUT,KMETCN,KDBANK, TIMES, NL, 

1 TREF , NDES , MAXNP , MNOPLY , PO , PI , DFH , NF , CODEH , MDES , TOLD , 

2 TNEW, MSH , NSH , QSH , FPI , MPI , DPI , FSIGI , MSIGI , DSIGI , FDOTI , 

3 MDOTI , DDOTI , FTVCI , MTVCI , DTVCI , TEMPER, NMECHC , NTHERC , 

4 PINTER, NLS) 

CALL HTMET ( A , MAXLEN , LD I TER , LD ALOW , TOLMET ) 

IF (EL3 .OR. EL7 .OR. EL11) THEN 
CALL MCPROP (IK, KOUT , KPOST 0 , PC ) 

ELSE IF (EL75) THEN 

CALL MCPROP ( IK, KOUT , KPOSTO , PC) 

CALL PLYH (IK, KOUT, NLEP,THHF, CODEH, DFH, MDES, MNOPLY, MAXNP, PLH, FPI, 

1 MPI , ANGLG ) 

CALL PLPROP (IK, KOUT, KPOST5,NL, MNOPLY, TEMP, PLH, THSTRN) 

CALL PROPH ( IK , KOUT , KPROP , IMAT , HTHI CK , NL , PLH , AS , CS , BS , ALPH , DEN ) 
END IF 
IREC = JREC 

CALL MPROP (IK, IREC, KOUT, KLAM,KPOST2,KPOST3,KPOST4,KPOST5,NL, 

1 MNOPLY, TNEW) 

JREC = JREC + NL 
30 CONTINUE 
KCODE = 1 


CALL GETIME (TT1, TT2 ,DTI ) 
DTI = DTI 

IF (FABRIC) GO TO 6000 


SUB. MHPREP SETS UP THE INPUT FILE KMHOST FOR THE MODULE SOLVEH 
(MHOST) . THE SUB. SOLVEH DETERMINES THE GLOBAL SOLUTION. THE 
STRESS RESULT. & DISP. ARE OBTAINED FROM THE OUTPUT FILES KMHFOR 
& KDISP. 

CALL MHPREP (SPARS, TEM, PRES, ANGULA, CFORCE, MODES, BUCK, EL3 , EL7 , ELIO, 

1 EL11,EL75, KOUT , KREAD , KMHOST , MESH , NETOT , NPTOT , MAXNP , 

2 NBOUND , GNP , KEI , GRIDP1 , GRIDP2 , ANVDIF, KAD , NCFNOD , 

3 NCFD IR , COFORC , TLAY , PD IFO , PDI FI , AS , CS , BS , ALPH , DEN , PC , 

4 NPRES S , PRE SUR , NTRAN S , TRANG , THSTRN ) 
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C 

C 

C 

C 

C 


CALL SOLVEH 
CALL MHOUTP 


( ISIZE , IWORK) 

(EL3 , EL7 , ELIO , EL11 , EL75 , KOUT, KMHFOR, KDISP , MAXNP , GNP , 
THHF , MSH,NSH, QSH , SIG Z , DI SVH , REACF , STRMOM , ETRAN ) 


CALL GETIME (TT1,TT2,DTI) 

DT2 = DTI 

CALL BAN1 (KOUT, ITER, INCR,NMC, IRAMP, 1,TIMEE) 


CJOE 

119 


121 


ITERATE OVER EACH NODE TO SET THE FLAGS FOR FIBER, MATRIX, & 
INTERPHASE PROP. FAILURE (SUB. MPROP) . ADD THE CURRENT STRESSES TO 
THE INITIAL STRESSES IN THE SUB. MSTRES. 

TOTTr* s 1 

DO 120 IK - 1, MAXNP 
PRINT *, ’NODE = * , IK 
NL = NLEP(IK) 

HTHICK * THHF (IK) 

PO = DABS (PDIFO( IK) ) 

PI = DABS ( PDIFI ( IK) ) 

DO 119 IL - 1 , NL 
TOLD (IL, IK) = TREF 
IF (TEM) THEN 

IF (EL3 .OR. EL 7 .OR. EL11) TOLD(IL,IK) = TPGNP(6, IK, 1) 

TOLD ( IL, IK) =0.0 

IF (.NOT. EL7 ) TOLD ( IL, IK) = THERMO ( IL, IK, 1) 

END IF 
CONTINUE 
NLS = 2 
CONTINUE 
CALL PMETCN 


( IK , IREC , JREC , INTER , KSTR , KOUT , KMETCN , KDBANK , TIMES , NL , 

1 TREF , NDES , MAXNP , MNOPLY , PO , PI , DFH , NF , CODEH , MDES , TOLD , 

2 TNEW,MSH,NSH,QSH, FPI ,MPI ,DPI , FSIGI , MSIGI ,DSIGI , FDOTI , 

3 MDOTI , DDOTI , FTVCI , MTVCI , DTVCI , TEMPER, NMECHC , NTHERC, 

4 PINTER, NLS) 

CALL HTMET ( A , MAXLEN , LD I TER , LD ALOW , TOLMET ) 

IREC = JREC 


IF (LDITER .LE. LDALOW ) WRITE (KOUT, 4010) IK, LDITER, NLS 
IF (LDITER .LE. LDALOW) GO TO 122 
WRITE (KOUT, 4000) IK, LDITER, NLS 
NLS = NLS + 1 

IF (NLS .GT. LSALOW ) GO TO 9000 
GO TO 121 
122 CONTINUE 

IF (EL3 .OR. EL7 .OR. EL11) THEN 
CALL MCPROP ( IK , KOUT , KPOSTO , PC ) 

ELSE IF (EL75 ) THEN 

CALL PLYH (IK, KOUT, NLEP, THHF, CODEH, DFH, MDES, MNOPLY, MAXNP, PLH, FPI, 
1 MPI , ANGLG) 

CALL PLPROP (IK, KOUT, KPOST5,NL, MNOPLY, TEMP, PLH, THSTRN) 

CALL PROPH (IK, KOUT, KPROP,IMAT, HTHICK, NL, PLH, AS, CS,BS,ALPH, DEN) 

END IF 
IREC = JREC 

CALL MPROP (IK, IREC, KOUT, KLAM,KPOST2,KPOST3,KPOST4,KPOST5,NL, 

1 MNOPLY, TNEW) 

IREC = JREC 
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CALL MSTRES ( IK , IREC , KOUT , KSTR , KS IG , KLAM , KPOST6 , KPOST7 , KPOST8 , 
1 KPOST9 , NL) 

JREC « JREC + NL 
120 CONTINUE 


C TOTAL DISPLACEMENTS AND PLY STRESSES AFTER THE FIRST LOAD STEP 

C 

DO 130 IK = 1, MAXNP 
DO 130 J = 1 , NELNFR 

DISP(J,IK) = DISP( J, IK) + DISVH( J , IK) 

130 CONTINUE 

CALL PSTRES (KOUT, KSIG,EL3,EL7, ELIO, EL11,EL75,NLEP, ITER, INCR, 

1 NMC, IRAMP, 1, l,TIMEE,DISP,TOTREA,MAXNP) 

C 

C PERFORM MODAL ANALYSIS 

C 

IF (MMODE .GT. 0) THEN 

CALL MODESH (IWORK, ISIZE, SPARS, TEM, PRES, ANGULA, CFORCE, MODES, BUCK, 

1 EL3 , EL7 , ELIO , ELI 1 , EL7 5 , KOUT , KREAD , KMHOST , KEIGEN , 

2 TIMEE, INCR, MESH, NETOT,NPTOT,MAXNP,GNP,KEI,NBOUND, 

3 GRIDP1 , GRIDP2 , ANVDIF , KAD , NCFNOD , NCFDIR , COFORC , TLAY , 

4 PDIFO, PDIFI , AS , CS , BS , ALPH, DEN , EIGNVA, EIGNVE , PC, 

5 NPRESS , PRESUR, NTRANS , TRANG , THSTRN ) 

END IF 


C 

C PERFORM BUCKLING ANALYSIS 

C 

IF (MBUCK .GT. 0) THEN 

CALL BUCKSH ( IWORK, ISIZE, SPARS, TEM, PRES, ANGULA, CFORCE, MODES, BUCK, 
1 EL3,EL7, ELIO, EL11,EL75, KOUT, KREAD, KMHOST, KEIGEN, 

1 TIMEE , INCR , ME SH , NETOT , NPTOT , MAXNP , GNP , KEI , NBOUND , 

2 GRIDP1,GRIDP2, ANVDIF, KAD, NCFNOD, NCFDIR, COFORC, TLAY, 

3 PDIFO, PDIFI , AS , CS , BS , ALPH, DEN, EIGNVA, EIGNVE , PC, 

4 NPRESS, PRESUR, NTRANS, TRANG, THSTRN) 

END IF 

C 

C GENERATE A PATRAN NEUTRAL FILE 

C 

CALL WPOST (KNODAL,KTRAN, MAXNP, DISP , NSH, MSH, QSH, TIMEE) 

CALL GETIME (TT1,TT2,DTI) 

DT3 = DTI + DT2 + DTI 
WRITE (KOUT, 3000) DTI, DT2 , DT3 
6000 CONTINUE 

IF (.NOT. RESTAR) KCODE = 1 
IF (.NOT. RESTAR) ACCTIM = 0.0 
MST = 0 
C 

C LOOP OVER EACH MECHANICAL & THERMAL CYCLE 

C 

CJOE DO 1000 NMC = 1,NMECHC 
DO 1000 NMC * 1,1 
NORAMP = NTISTP - 1 
IF (.NOT. RESTAR) MRAMP = 1 

IF (RESTAR .AND. INCR .LT. LINC) MRAMP = IRAMP 
IF (RESTAR .AND. INCR .EQ. LINC) MRAMP = IRAMP + 1 
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LOOP OVER NUMBER OP LOADING RAMPS (NORAMP) 

DO 200 IRAMP = MRAMP, NORAMP 
NSTPI = IRAMP 
NSTPI1 * IRAMP + 1 
TIMEIN ■ TISTPS (NSTPI) 

TIMEI1 * TISTPS ( NSTPI 1) 

DTIME = (TIMEI1 - TIMEIN) /FLOAT (LINC) 

CALL LOADIN ( PRES, TEM , ANGULA , CFORCE , LINC, MAXNP, NSTPI 1, NSTPI, NCFOR, 

1 NPRES , CFVAL , COFORC , TPGNP , PDIFO, PDIFI , 0 . 0 , TLAY , ANGVEL , 

2 ANVDIF,NTISTP, PREVAL,PRESUR) 


C INCREMENT OVER EACH RAMP ACCORDING TO NO. OF LOAD INCR. (LINC) 

C 

IF (.NOT. RESTAR) MINCR = 1 

IF (RESTAR) MINCR = INCR + 1 

IF (RESTAR .AND. MINCR .GT. LINC) MINCR = 1 

DO 500 INCR = MINCR, LINC 

RLINC = FLOAT (LINC) 

RINCR = FLOAT (INCR) 

MST = MST + 1 

CALL GETIME (TT1,TT2,DTI) 

DT7 = DTI 
TTIME1 — TT1 

TIMEE = TIMEIN + RINCR*DTIME 
TIMES (1) » TIMEIN + FLOAT ( INCR- 1 ) *DTIME 
TIMES (2) = TIMEE 
C 

C ITERATE OVER EACH LOAD INCREMENT (INCR) 

C 

ITER = 0 
600 CONTINUE 
DTIM = 0 

CALL GETIME (TT1,TT2,DTI) 

ITER = ITER + 1 


C IF NO. OF ITERATIONS ( ITER ) EXCEEDS MITER WRITE A RESTART FILE 

C 

IF (ITER .LE. MITER) GO TO 7000 
WRITE (KOUT, 3100 ) IRAMP, INCR 

IF (INCR .EQ. 1 .AND. IRAMP .EQ. 1) GO TO 7100 
IF (INCR .EQ. 1) IINCR = LINC 
IF (INCR .GT. 1) IINCR = INCR - 1 
I IRAMP = IRAMP 

IF (INCR .EQ. 1 .AND. IRAMP .GT. 1) IIRAMP = IRAMP - 1 
WRITE (KOUT, 3200) IIRAMP, IINCR 
GO TO 9000 
7100 CONTINUE 

WRITE (KOUT, 3300) 

GO TO 1000 
7000 CONTINUE 


IF ( KCODE .LT. 1) GO TO 7050 

CALL MHPREP (SPARS, TEM, PRES, ANGULA, CFORCE, MODES, BUCK, EL3,EL7, ELIO, 
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1 

2 

3 

4 

CALL SOLVEH 
CALL MHOUTP 


1 

CALL GETIME 
DT9 = DTI 
7050 CONTINUE 


EL11, EL7 5 , KOUT , KRE AD , KMHOST , MESH , NETOT , NPTOT , MAXNP , 

N BOUND , GNP , KEI , GRIDP1 , GRIDP2 , ANVDIF , KAD , NCFNOD , 
NCFDIR , COFORC , TLAY , PD I FO , PD I FI , AS , CS , BS , ALPH , DEN , PC , 
NPRESS , PRESUR, NTRANS , TRANG , THSTRN ) 

( I S I ZE / I WORK ) 

( EL3 , EL 7 , ELIO , EL11 , EL7 5 , KOUT , KMHFOR, KDISP , MAXNP , GNP , 
THHF , MSH , NSH , QSH , SIGZ / DI SVH , RE ACF , STRMOM , ETRAN ) 

( TT1 , TT2 , DTI ) 


Tprr = 1 

DO 245 IK = 1, MAXNP 
HTHICK = THHF (IK) 

NL = NLEP(IK) 

DO 265 IL = 1,NL 
TOLD (IL, IK) = TREF 
TNEW(IL, IK) = TREF 
TEMP (IL, IK) =0.0 
IF (TEM) THEN 

IF (EL3 .OR. EL 7 .OR. EL11) THEN 

DTEMP = (TPGNP(6, IK,NSTPI1) - TPGNP (6, IK,NSTPI) ) /RLINC 
TOLD ( IL, IK) = TPGNP ( 6 , IK, NSTPI ) + (RINCR - l.)*DTEMP 
TNEW( IL, IK) = TPGNP (6, IK, NSTPI) + RINCR*DTEMP 

TEMP( IL, IK) = (TPGNP (6, IK, NSTPI) - TPGNP ( 6, IK, NSTPI1 )) /RLINC 
ELSE 

DTEMP = ( THERMO (IL, IK, NSTPI1) - THERMO( IL, IK, NSTPI )) /RLINC 
TOLD ( IL, IK) = THERMO (IL, IK, NSTPI) + (RINCR - l.)*DTEMP 
TNEW(IL, IK) = THERMO (IL, IK, NSTPI) + RINCR*DTEMP 

TEMP(IL,IK) = ( THERMO (IL, IK, NSTPI) - THERMO( IL, IK, NSTPI1) ) /RLINC 
END IF 
END IF 
265 CONTINUE 

PO = DABS(PDIFO(IK) ) 

PI = DABS ( PDIFI ( IK) ) 

NLS = 2 
221 CONTINUE 


CALL PMETCN ( IK, IREC, JREC, INTER, KSTR, KOUT, KMETCN,KDBANK, TIMES, NL, 

1 TREF , NDES , MAXNP , MNOPLY , PO , PI , DFH , NF, CODEH , MDES , TOLD , 

2 TNEW , MSH , NSH , QSH , FPI , MPI , DPI , FSIGI , MSIGI , DSIGI , FDOTI , 

3 MDOTI , DDOTI , FTVCI , MTVCI , DTVCI , TEMPER, NMECHC , NTHERC , 

4 PINTER, NLS) 

CALL HTMET ( A , MAXLEN , LD ITER, LD ALOW , TOLMET ) 

IREC = JREC 

C IF (LDITER .LE. LDALOW ) WRITE (KOUT, 4010) IK, LDITER, NLS 

IF (LDITER .LE. LDALOW) GO TO 222 
WRITE (KOUT, 4000) IK, LDITER, NLS 
NLS = NLS + 1 

IF (NLS .GT. LSALOW ) GO TO 9000 
GO TO 221 
222 CONTINUE 

CALL GETIME (TT1,TT2,DTI) 

DTIM = DTIM + DTI 

IF (EL3 .OR. EL7 .OR. EL11) THEN 

CALL MCPROP ( IK , KOUT , KPOSTO , PC ) 

ELSE IF (EL75 ) THEN 
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CALL PLYH ( IK , KOUT , NLEP , THHF , CODEH , DFH , MDES , MNOPLY , MAXNP , PLH, FPI , 
1 MPI, ANGLG) 

CALL PLPROP (IK, KOUT , KPOST 5 , NL , MNOPLY , TEMP , PLH , THSTRN ) 

CALL PROPH ( IK , KOUT , KPROP , IMAT , HTHICK , NL , PLH , AS , CS , BS , ALPH , DEN ) 
END IF 

JREC = JREC + NL 
245 CONTINUE 

CALL PSTRES (KOUT, KSIG, EL3 , EL7 , ELIO , ELI 1 , EL75 , NLEP , ITER, INCR, 

1 NMC,IRAMP,5,5,TIMEE,DISVH,REACF,MAXNP) 


C 

C CHECK FOR CONVERGENCE 

C 

CALL CHECKD (CONV, LINEAR, KOUT, KCODE, TOL, NMC, IRAMP, INCR, ITER, NPTOT, 
1 D I S VH , DERORO , NELNFR ) 

KCODE = 1 

CALL GETIME (TT1,TT2,DTI) 

DT11 = DT9 + DT10 + DTIM 
WRITE (KOUT, 2007) DT9 ,DTIM,DT11 
IF (CONV) GO TO 280 
GO TO 600 
280 CONTINUE 


C 

C 

C 


IF CONVERGED, DETER. NEW FIBER, MATRIX, INTERPHASE, AND STRENGTHS 
JREC — 1 

DO 400 IK = 1, MAXNP 
HTHICK = THHF (IK) 

NL = NLEP (IK) 

DO 365 IL = 1,NL 
TOLD (IL, IK) = TREF 
TNEW ( IL, IK) = TREF 
TEMP (IL, IK) =0.0 
IF (TEM) THEN 

IF (EL3 .OR. EL7 .OR. EL11) THEN 

DTEMP = (TPGNP ( 6, IK, NSTPI1 ) - TPGNP (6, IK,NSTPI ) ) /RLINC 

TOLD ( IL, IK) = TPGNP ( 6 , IK, NSTPI ) + (RINCR - l.)*DTEMP 

TNEW( IL, IK) = TPGNP ( 6 , IK, NSTPI ) + RINCR*DTEMP 

TEMP(IL, IK) = (TPGNP (6, IK, NSTPI) - TPGNP (6, IK, NSTPI1 )) /RLINC 

ELSE 

DTEMP = ( THERMO ( IL, IK, NSTPI 1 ) - THERMO ( IL, IK, NSTPI )) /RLINC 
TOLD (IL, IK) » THERMO (IL, IK, NSTPI) + (RINCR - l.)*DTEMP 
TNEW(IL, IK) = THERMO ( IL, IK, NSTPI ) + RINCR*DTEMP 

TEMP (IL, IK) = (THERMO (IL, IK, NSTPI) - THERMO ( IL, IK, NSTPI1 )) /RLINC 
END IF 
END IF 
365 CONTINUE 

PO = DABS ( PDIFO( IK) ) 

PI = DABS (PDIFI ( IK) ) 

NLS * 2 
421 CONTINUE 

CALL PMETCN 


1 

2 

3 

4 


( IK , IREC , JREC , INTER , KSTR, KOUT , KMETCN , KDBANK , TIMES , NL , 
TREF , NDES , MAXNP , MNOPLY , PO, PI , DFH , NF , CODEH , MDES , TOLD , 
TNEW, MSH,NSH,QSH, FPI, MPI, DPI, FSIGI , MSIGI , DSIGI , FDOTI , 
MDOTI , DDOTI , FTVCI , MTVCI , DTVCI , TEMPER, NMECHC, NTHERC , 
PINTER, NLS) 
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CALL HTMET ( A , MAXLEN , LD ITER , LD ALOW , TOLMET ) 

IREC = JREC 

C IF (LDITER .LE. LDALOW ) WRITE (KOUT,4010) IK, LDITER, NLS 

IF (LDITER .LE. LDALOW) GO TO 422 
WRITE (KOUT, 4000) IK, LDITER, NLS 
NLS = NLS + 1 

IF (NLS .GT. LSALOW ) GO TO 9000 

GO TO 421 
422 CONTINUE 

CALL MPROP (IK, IREC, KOUT, KLAM,KPOST2,KPOST3,KPOST4,KPOST5,NL, 

1 MNOPLY , TNEW ) 

IREC = JREC 

CALL MSTRES ( IK, IREC, KOUT, KSTR,KSIG,KLAM,KPOST6,KPOST7 , KPOST8, 

1 KPOST9 , NL) 

JREC = JREC + NL 
400 CONTINUE 
C 

C UPDATE THE TOTAL DISPL. (DISP) WITH INCR. DISPL. (DISVH) 

C 

DO 415 IK = 1,MAXNP 
DO 415 J = 1 , NELNFR 

DISP ( J, IK) = DISP( J, IK) + DISVH( J , IK) 

415 CONTINUE 

CALL PSTRES (KOUT, KSIG, EL3 , EL7 , ELIO, EL11 ,EL75,NLEP, ITER, INCR, 

1 NMC , IRAMP ,3,3, TIMEE,DISP, TOTRE A , MAXNP ) 

CALL GETIME (TT1,TT2,DTI) 

DTI = TT1 - TTIME1 
WRITE (KOUT, 2040) DTI 
ACCTIM = TT2 
C 

C PERFORM MODAL ANALYSIS 

C 

CALL MODESH (IWORK, ISIZE, SPARS, TEM, PRES, ANGULA,CFORCE, MODES, BUCK, 

1 EL3 , EL7 , ELIO , ELI 1 , EL75 , KOUT , KREAD , KMHOST , KEIGEN , 

2 TIMEE, INCR, MESH, NETOT,NPTOT, MAXNP, GNP,KEI,NBOUND, 

3 GRIDP1, GRIDP2 , ANVDIF , KAD , NCFNOD , NCFDIR, COFORC, TLAY , 

4 PDIFO,PDIFI , AS, CS,BS, ALPH,DEN,EIGNVA,EIGNVE, PC, 

5 NPRESS , PRESUR , NTRANS , TRANG , THSTRN ) 

C 

C PERFORM BUCKLING ANALYSIS 

C 

CALL BUCKSH ( IWORK, ISIZE, SPARS, TEM, PRES, ANGULA, CFORCE, MODES, BUCK, 
1 EL3 , EL7 , ELIO , ELI 1 , EL7 5 , KOUT , KREAD , KMHOST , KEIGEN , 

1 TIMEE , INCR , MESH , NETOT , NPTOT , MAXNP , GNP , KEI , NBOUND , 

2 GRIDP1 , GRIDP2 , ANVDIF, KAD , NCFNOD , NCFDIR, COFORC, TLAY , 

3 PDIFO, PDIFI , AS , CS , BS , ALPH , DEN, EIGNVA, EIGNVE , PC , 

4 NPRES S , PRE SUR , NTRANS , TRANG , THSTRN ) 

C 

C GENERATE A PATRAN NEUTRAL FILE 

C 

CALL WPOST (KNODAL,KTRAN, MAXNP, DISP, NSH,MSH,QSH, TIMEE) 

C 

IF (MST .LT. MSTART) GO TO 500 
I IRAMP = IRAMP 
I INCR = INCR 
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C 

C WRITE DATA TO RESTART FILE KREAD 

C 

9000 CONTINUE 
JCODE = 1 

ACCTIM * ACCTIM + ACCT 

CALL RESTART ( JCODE, KREAD, KSTART, KLAM,KSIG,KSTR, NMC, NORAMP, I RAMP, 
1 INCR, ACCTIM, NLEP , DISP , TOTREA, MAXNP , NELNFR ) 

GO TO 1000 


C 

JCODE = 1 

ACCTIM = ACCTIM + ACCT 

CALL RESTART (JCODE, KREAD, KSTART, KLAM,KSIG,KSTR, NMC, NORAMP, IRAMP, 
1 INCR , ACCTIM , NLEP , D I SP , TOTREA, MAXNP , NELNFR ) 

500 CONTINUE 
200 CONTINUE 

C C ALL PSTRES (KOUT,KSIG,EL3 , EL7, ELIO, EL11 ,EL75 , NLEP, ITER, INCR, 

C 1 NMC, IRAMP, 3, 4, TIMEE, DISP, TOTREA, MAXNP) 

1000 CONTINUE 

CALL GETIME (TTl,TT2,DTI) 

JtCCTIM = ACCTIM + DTI 
WRITE (KOUT, 2150) ACCTIM 


C 

2007 FORMAT 
1 
2 
3 

2018 FORMAT 
1 
2 
3 

2040 FORMAT 
1 

2150 FORMAT 
1 

2160 FORMAT 
1 
2 
3 

3000 FORMAT 
1 
1 
2 

3100 FORMAT 
1 

3200 FORMAT 
1 

3300 FORMAT 

4000 FORMAT 
1 

4010 FORMAT 
1 


( / 1 TIME REQUIRED TO : CALCULATE THE STRESS ’ , 

’RESULTANTS IN MHOST ’ , 6X, F10 . 3 , ’ SEC. ’,/, 20X, ’ DETERMINE ’ , 
’ PLY AND NODAL PROPERTIES ’,6X,F10.3,’ SEC.’,/,20X, 

’PERFORM THIS ITERATION ’ ,14X,F10.3, ’ SEC.’) 

(/’ TIME REQUIRED TO : READ IN DATA FROM ’, 

’RESTART FILE ’,F10.3,’ SEC.’,/, 

2 OX, ’TOTAL TIME REQUIRED FOR THE ANALYSIS SO FAR’, 

3X, F10. 3 , ’ SEC.’) 

(/’ TIME REQUIRED TO EVALUATE THIS LOAD INCREMENT ’, 20X, 
F10. 3 , ’ SEC.’) 

(/’ TIME REQUIRED TO CARRY OUT THE ANALYSIS 
F10.3, ’ SEC.’) 

( / ’ TIME REQUIRED TO : ’ , 

’PERFORM FIRST LOAD STEP ’ , 17X,F10.3, ’ SEC. 

’DO THE NONLINEAR ANALYSIS ’ , 17X,F10.3, ’ SEC. 

’CARRY OUT THE ENTIRE ANALYSIS' , 17X,F10. 3, ’ SEC. 

(/’ TIME REQUIRED TO : ' , 

'LOOP THRU METCAN ' , 17X, F10 . 3, ’ SEC. 

'DETER. THE DISPL. IN MHOST ' , 17X, F10 . 3, ’ SEC. 

'ANALYZE THE FIRST LOAD STEP ’ , 17X, F10. 3 , ’ SEC.') 

(/' THE ANALYSIS HAS FAILED TO CONVERGE FOR RAMP NO.', IX, 
13,! AND LOAD INCREMENT NO.’, IX, 13) 

(/' A RESTART FILE HAS BEEN CREATED FOR RAMP NO.',5X,I3, 

' AND LOAD INCREMENT NO.', IX, 13) 

( / ' A RESTART FILE WAS NOT CREATED ' ) 

(' NUMBER OF ITERATIONS FOR NODE # ',14,' IS ’,12,' WITH 
12,' LOAD STEPS THAT METCAN FAILED TO CONVERGE ') 

(' NUMBER OF ITERATIONS REQUIRED FOR NODE # ',14,' IS ', 
12,' WITH ',12,' LOAD STEPS TO HAVE CONVERGENCE ') 


, 20X, 


' , / , 20X, 
1 ,/,20X, 
' ) 

' ,/,20X, 
' ,/,20X, 


RETURN 

END 
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